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Abstract 

ISAJET is a Monte Carlo program which simulates pp, pp and e + e~ interactions at high 
energies. This document summarizes the physics underlying the program and describes how 
to use it. New features of Version 7.40 include generation of W + H and Z + H events, 
non-minimal GMSB models, and use of the experimental a s (Mz) by default. 
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1 Introduction 



ISAJET is a Monte Carlo program which simulates pp and pp interactions at high energies. 
It also generates e + e~ scattering events, although this is less developed. ISAJET is based 
on perturbative QCD plus phenomenological models for parton and beam jet fragmentation. 
Events are generated in four distinct steps: 

• A primary hard scattering is generated according to the appropriate QCD cross section. 

• QCD radiative corrections are added for both the initial and the final state. 

• Partons are fragmented into hadrons independently, and particles with lifetimes less 
than about 1CT 12 seconds are decayed. 

• Beam jets are added assuming that these are identical to a minimum bias event at the 
remaining energy. 

ISAJET incorporates ISASUSY, which evaluates branching ratios for the minimal su- 
persymmetric extension of the standard model. H. Baer and X. Tata are coauthors of this 
package, and they have done the original calculations with various collaborators. See the 
ISASUSY documentation in the patch Section [T^. 

ISAJET is supported for ANSI Fortran and for Cray, DEC Ultrix, DEC VMS, HP/9000 
7xx, IBM VM/CMS 370 and 30xx, IBM AIX RS/6000, Linux, Silicon Graphics 4D, and Sun 
computers. The CDC 7600 and ETA 10 versions are obsolete and are no longer supported. 
It is written mainly in ANSI standard FORTRAN 77, but it does contain some extensions 
except in the ANSI version. The code is maintained with a combination of RCS, the Revision 
Control System, and the Patchy code management system, which is part of the CERN 
Library. The original sources are kept on physgi01.phy.bnl.gov in "isajet/ isalibrary/RCS; 
decks revised in release n.nn are kept in ~isajet/isalibrary/nnn. ISAJET is supplied to 
BNL, CERN, Fermilab, and SLAC; it is also available by anonymous ftp from 

\protect\vrule widthOpt\protect\href {f tp : //penguin . phy . bnl . gov/pub/ isaj et}{f tp : //penguin 

or by request from the authors. 

Patch ISAPLT contains the skeleton of an HBOOK histogramming job, a trivial calorime- 
ter simulation, and a jet-finding algorithm. (The default is HBOOK4; HBOOK3 can be 
selected with a Patchy switch.) These are provided for convenience only and are not sup- 
ported. 
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2 Physics 



ISAJET is a Monte Carlo program which simulates pp and pp interactions at high energy. 
(It also contains an e + e~ event generator, but this is much less developed.) The program 
incorporates perturbative QCD cross sections, initial state and final state QCD radiative 
corrections in the leading log approximation, independent fragmentation of quarks and gluons 
into hadrons, and a phenomenological model tuned to minimum bias and hard scattering 
data for the beam jets. 

2.1 Hard Scattering 

The first step in simulating an event is to generate a primary hard scattering according to 
some QCD cross section. This has the general form 

a = a F(x 1 ,Q 2 )F(x 2 ,Q 2 ) 

where o"o is a cross section calculated in QCD perturbation theory, F(x, Q 2 ) is a struc- 
ture function incorporating QCD scaling violations, X\ and x 2 are the usual parton model 
momentum fractions, and Q 2 is an appropriate momentum transfer scale. 

For each of the processes included in ISAJET, the basic cross section a is a two-body 
one, and the user can set limits on the kinematic variables and type for each of the two 
primary jets. For DRELLYAN and WPAIR events the full matrix element for the decay of 
the W's into leptons or quarks is also included. 

The following processes are available: 

2.1.1 Minbias 

No hard scattering at all, so that the event consists only of beam jets. Note that at high 
energy the jet cross sections become large. To represent the total cross section it is better 
to use a sample of TWOJET events with the lower limit on pt chosen to give a cross section 
equal to the inelastic cross section or to use a mixture of MINBIAS and TWOJET events. 

2.1.2 Twojet 

All order a 2 QCD processes, which give rise in lowest order to two high-p t jets. Included 
are, e.g. 

9 + 9 -> 9 + 9 

g + q -»• g + q 
g + g -»• q + q 

Masses are neglected for c and lighter quarks but are taken into account for b and t quarks. 
The Q 2 scale is taken to be 

Q 2 = 2stu/{s 2 + t 2 + u 2 ) 

The default parton distributions are those of the CTEQ Collaboration, fit CTEQ3L, using 
lowest order QCD evolution. Two older fits, Eichten, Hinchliffe, Lane and Quigg (EHLQ), 
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Set 1, and Duke and Owens, Set 1, are also included. There is also an interface to the CERN 
PDFLIB compilation of parton distributions. Note that structure functions for heavy quarks 
are included, so that processes like 

g + t^g + t 

can be generated. The Duke-Owens parton distributions do not contain b or t quarks. 

Since the t is so heavy, it decays before it can hadronize, so instead of t hadrons a t 
quark appears in the particle list. It is decayed using the V — A matrix element including 
the W propagator with a nonzero width, so the same decays should be used for m t < mw 
and m t > mw', the W should not be listed as part of the decay mode. The partons are then 
evolved and fragmented as usual; see below. The real or virtual W and the final partons 
from the decay, including any radiated gluons, are listed in the particle table, followed by 
their fragmentation products. Note that for semileptonic decays the leptons appear twice: 
the lepton parton decays into a single particle of the same type but in general somewhat 
different momentum. In all cases only particles with IDCAY = should be included in the 
final state. 

A fourth generation x, y is also allowed. Fourth generation quarks are produced only 
by gluon fusion. Decay modes are not included in the decay table; for a sequential fourth 
generation they would be very similar to the t decays. In decays involving quarks, it is 
essential that the quarks appear last. 

2.1.3 Drellyan 

Production of a W in the standard model, including a virtual 7, a W + , a W~, or a Z°, and 
its decay into quarks or leptons. If the transverse momentum QTW of the W is fixed equal 
to zero then the process simulated is 

q + q^W -> q + q 
-> l + i 

Thus the W has zero transverse momentum until initial state QCD corrections are taken 
into account. If non-zero limits on the transverse momentum q t for the W are set, then 
instead the processes 

q + q -> W + g 
g + q -> W + q 

are simulated, including the full matrix element for the W decay. These are the dominant 
processes at high q t , but they are of course singular at q t — 0. A cutoff of the 1/qf singularity 
is made by the replacement 

!/«?-> l/\fof+qJo g t 2 = (.2GeV)M 

This cutoff is chosen to reproduce approximately the qt dependence calculated by the sum- 
mation of soft gluons and to give about the right integrated cross section. Thus this option 
can be used for low as well as high transverse momenta. 
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The scale for QCD evolution is taken to be proportional to the mass for lowest order 
Drell-Yan and to the transverse momentum for high-p t Drell-Yan. The constant is adjusted 
to get reasonable agreement with the W + n jet cross sections calculated from the full QCD 
matrix elements by F.A. Berends, et al., Phys. Lett. B224, 237 (1989). 

For the processes g + b — > W + t and g + t — > Z + t, cross sections with a non-zero top 
mass are used for the production and the W/Z decay. These were calculated using FORM 
1.1 by J. Vermaseren. The process g + t ^ W + b is not included. Both g + b — > W~ + t and 
g + t — > W~ + 6 of course give the same W~ + t + b final state after QCD evolution. While 
the latter process is needed to describe the m t = 0(!) mass singularity for q t ^> m t , it has a 
pole in the physical region at low qt from on-shell t W + b decays. There is no obvious 
way to avoid this without introducing an arbitrary cutoff. Hence, selecting only W + b will 
produce a zero cross section. The Q 2 scale for the parton distributions in these processes 
is replaced by Q 2 + m 2 ; this seems physically sensible and prevents the cross sections from 
vanishing at small q t . 

2.1.4 Photon 

Single and double photon production through the lowest order QCD processes 

g + q -> 7 + <? 
q + q -> 7 + 3 
q + q -> 7 + 7 

Higher order corrections are not included. But 7's, W's, and Z's are radiated from final 
state quarks in all processes, allowing study of the bremsstrahlung contributions. 

2.1.5 Wpair 

Production of pairs of W bosons in the standard model through quark-antiquark annihilation, 

q + q -> W + + W~ 
-> Z° + Z° 

-»• W^ + + 7,H/-+7 

The full matrix element for the W decays, calculated in the narrow resonance approximation, 
is included. However, the higher order processes, e.g. 

q + q -> q + q + W + + W~ 

are ignored, although they in fact dominate at high enough mass. Specific decay modes can 
be selected using the WMODEi keywords. 
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2.1.6 Higgs 

Production and decay of the standard model Higgs boson. The production processes are 

g + g — > H (through a quark loop) 

q + q — > H (with t + i dominant) 

W + + W~ — > if (with longitudinally polarized W) 

Z° + Z° -> if (with longitudinally polarized Z) 

If the (Standard Model) Higgs is lighter than 2M W) then it will decay into pairs of fermions 
with branching ratios proportional to mj. If it is heavier than 2Mw, then it will decay 
primarily into W + W~ and Z°Z° pairs with widths given approximately by 



r(H -> z°z°) 

Numerically these give approximately 

= 0.5 TeV 



G F Mfj 
8tt^2 
G F Mfj 



M H ^ 



3 



ITeV 



The width proportional to arises from decays into longitudinal gauge bosons, which like 
Higgs bosons have couplings proportional to mass. 

Since a heavy Higgs is wide, the narrow resonance approximation is not valid. To obtain 
a cross section with good high energy behavior, it is necessary to include a complete gauge- 
invariant set of graphs for the processes 



w + w~ 


-> w + w 


w + w~ 




z°z° 




z°z° 





with longitudinally polarized H^~, and Z° bosons in the initial state. This set of 

graphs and the corresponding angular distributions for the W + , W~, and Z° decays have 
been calculated in the effective W approximation and included in HIGGS. The W structure 
functions are obtained by integrating the EHLQ parameterization of the quark ones term by 
term. The Cabibbo-allowed branchings 

q -> W + + q' 
q -> W~ + q' 
q -> Z° + q 

are generated by backwards evolution, and the standard QCD evolution is performed. This 
correctly describes the W collinear singularity and so contains the same physics as the 
effective W approximation. 
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If the Higgs is lighter than 2M W , then its decay to 77 through W and t loops may be 
important. This is also included in the HIGGS process and may be selected by choosing GM 
as the jet type for the decay. 

If the Higgs has Mz < Mh < 2Mz, then decays into one real and one virtual Z° are 
generated if the ZO ZO decay mode is selected, using the calculation of Keung and Marciano, 
Phys. Rev. D30, 248 (1984). Since the calculation assumes that one Z° is exactly on shell, 
it is not reliable within of order the Z° width of M H = 2M Z ; Higgs and and Z°Z° masses 
in this region should be avoided. The analogous Higgs decays into one real and one virtual 
charged W are not included. 

Note that while HIGGS contains the dominant graphs for Higgs production and graphs 
for W pair production related by gauge invariance, it does not contain the processes 

q + q -> W + W~ 
q + q -> Z°Z° 

which give primarily transverse gauge bosons. These must be generated with WPAIR. 

If the MSSMi or SUGRA keywords are used with HIGGS, then one of the three MSSM neutral 
Higgs is generated instead using gluon-gluon and quark-antiquark fusion with the appropriate 
SUSY couplings. Since heavy CP even SUSY Higgs are weakly coupled to W pairs and CP 
odd ones are completely decoupled, WW fusion and WW — > WW scattering are not included 
in the SUSY case. (WW — > WW can be generated using the Standard Model process with 
a light Higgs mass, say 100 GeV.) The MSSM Higgs decays into both Standard Model and 
SUSY modes as calculated by ISASUSY are included. For more discussion see the SUSY 
subsection below and the writeup for ISASUSY. The user must select which Higgs to generate 
using HTYPE; see Section 6 below. If a mass range is not specified, then the range mass 
M H ± hY H is used by default. (This cannot be done for the Standard Model Higgs because 
it is so wide for large masses.) Decay modes may be selected in the usual way. 

2.1.7 WHiggs 

Generates associated production of gauge and Higgs bosons, i.e., 

q + q^H + W,H + Z, 

in the narrow resonance approximation. The desired subprocesses can be selected with 
JETTYPEi, and specific decay modes of the W and/or Z can be selected using the WMODEi 
keywords. Standard Model couplings are assumed unless SUSY parameters are specified, in 
which case the SUSY couplings are used. 

2.1.8 SUSY 

Generates pairs of supersymmetric particles from gluon-quark or quark-antiquark fusion. If 
the MSSMi or SUGRA parameters defined in Section 6 below are not specified, then only 
gluinos and squarks are generated: 

9 + 9 -> 9 + 9 
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g + g -> g + g 
g + q -> # + g 

# + # -> <? + g~ 
g + g — ► g + g 
g + g — > g + g 

Left and right squarks are distinguished but assumed to be degenerate. Masses can be 
specified using the GAUGINO, SQUARK, and SLEPTON parameters described in Section 6. No 
decay modes are specified, since these depend strongly on the masses. The user can either 
add new modes to the decay table (see Section 9) or use the FORCE or F0RCE1 commands 
(see Section 6). 

If MSSMA, MSSMB, and MSSMC are specified, then the ISASUSY package is used to calculate 
the masses and decay modes in the minimal supersymmetric extension of the standard model 
(MSSM), assuming SUSY grand unification constraints in the neutralino and chargino mass 
matrix but allowing some additional flexibility in the masses. The scalar particle soft masses 
are input via MSSMi, so that the physical masses will be somewhat different due to D-term 
contributions and mixings for 3rd generation sparticles. t\ and ti production and decays are 
now included. The lightest SUSY particle is assumed to be the lightest neutralino Z\. If the 
MSSMi parameters are specified, then the following additional processes are included using 
the MSSM couplings for the production cross sections: 

g + q -> Zi + q, W { + q 

q + q -> Zi + g, W t + g 

q + q -> Wi + Zj 

q + q -> W+ + W- 

q + q -> Zi + Zj 

q + q -> £ + + £~, v + v 

Processes can be selected using the optional parameters described in Section 6 below. 

An optional keyword MSSMD can be used to specify the second generation masses, which 
otherwise are assumed degenerate with the first generation. An optional keyword MSSME can 
be used to specify values of the U(l) and 577(2) gaugino masses at the weak scale rather 
than using the default grand unification values. The chargino and neutralino masses and 
mixings are then computed using these values. 

Instead of using the MSSMi parameters, one can use the SUGRA parameter to specify in the 
minimal supergravity framework. This assumes that the gauge couplings unify at a GUT 
scale and that SUSY breaking occurs at that scale with universal soft breaking terms, which 
are related to the weak scale using the renormalization group. The parameters of the model 
are 

• m : the common scalar mass at the GUT scale; 

• mi/2- the common gaugino mass at the GUT scale; 
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• A : the common soft trilinear SUSY breaking parameter at the GUT scale; 

• tan/5: the ratio of Higgs vacuum expectation values at the electroweak scale; 

• sgn/i = ±1: the sign of the Higgsino mass term. 

The renormalization group equations are solved iteratively to determine all the electroweak 
SUSY parameters from these data assuming radiative electroweak symmetry breaking but 
not other possible constraints such as b-tau unification or limits on proton decay. 

The assumption of universality at the GUT scale is rather restrictive and may not be 
valid. A variety of non-universal SUGRA (NUSUGRA) models can be generated using the 
NUSUG1, . . . , NUSUG5 keywords. These might be used to study how well one could test the 
minimal SUGRA model. 

An alternative to the SUGRA model is the Gauge Mediated SUSY Breaking (GMSB) 
model of Dine, Nelson, and collaborators. In this model SUSY breaking is communicated 
through gauge interactions with messenger fields at a scale M m small compared to the Planck 
scale and are proportional to gauge couplings times A m . The messenger fields should form 
complete SU(5) representations to preserve the unification of the coupling constants. The 
parameters of the GMSB model, which are specified by the GMSB keyword, are 

• A m = F m /M m : the scale of SUSY breaking, typically 10-100 TeV; 

• M m > A m : the messenger mass scale; 

• N 5 : the equivalent number of 5 + 5 messenger fields. 

• tan/5: the ratio of Higgs vacuum expectation values at the electroweak scale; 

• sgn/x = ±1: the sign of the Higgsino mass term; 

• Cgrav > 1: the ratio of the gravitino mass to the value it would have had if the only 
SUSY breaking scale were F m . 

In GMSB models the lightest SUSY particle is always the nearly massless gravitino G. The 
parameter C grav scales the gravitino mass and hence the lifetime of the next lightest SUSY 
particle to decay into it. The NOGRAV keyword can be used to turn off gravitino decays. 

A variety of non-minimal GMSB models can be generated using additional parameters 
set with the GMSB2 keyword. These additional parameters are 

• ft, an extra factor multiplying the gaugino masses at the messenger scale. (Models 
with multiple spurions generally have ft < 1.) 

• 5Mfj d , 5Mfj u , Higgs mass-squared shifts relative to the minimal model at the messenger 
scale. (These might be expected in models which generate \x realistically.) 

• D Y (M), a U(l)y messenger scale mass-squared term (D-term) proportional to the 
hyper charge Y . 
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• N 5l , N$ 2 , and N 53 , independent numbers of gauge group messengers. They can be 
non-integer in general. 

For discussions of these additional parameters, see S. Dimopoulos, S. Thomas, and J.D. 
Wells, |hep-ph/ 960933% Nucl. Phys. B488, 39 (1997), and S.P. Martin, |hep-ph/ 9608223 , 
Phys. Rev. D55, 3177 (1997). 

Gravitino decays can be included in the general MSSM framework by specifying a grav- 
itino mass with MGVTNO. The default is that such decays do not occur. 

The ISASUSY program can also be used independently of the rest of ISAJET, either 
to produce a listing of decays or in conjunction with another event generator. Its physics 



assumptions are described in more detail in Section [12]. The ISASUGRA program can also 
be used independently to solve the renormalization group equations with SUGRA, GMSB, 
or NUSUGRA boundary conditions and then to call ISASUSY to calculate the decay modes. 

Generally the MSSM, SUGRA, or GMSB option should be used to study supersymmetry 
signatures; the SUGRA or GMSB parameter space is clearly more manageable. The more 
general option may be useful to study alternative SUSY models. It can also be used, e.g., to 
generate pointlike color-3 leptoquarks in technicolor models by selecting squark production 
and setting the gluino mass to be very large. The MSSM or SUGRA option may also be 
used with top pair production to simulate top decays to SUSY particles. 

2.1.9 e+e" 

An e + e~ event generator is included in ISAJET, although it is less well developed than some 
others. The Standard Model processes included are e + e~ annihilation through 7 and Z to 
quarks and leptons, and production of W + W~ and Z°Z° pairs. In contrast to WPAIR and 
HIGGS for the hadronic processes, the produced W's and Z's are treated as particles, so 
their spins are not properly taken into account in their decays. (Because the PU's and Z's 
are treated as particles, their decay modes must be selected using FORCE or F0RCE1, not 
WMODEi. See Section [6] below.) Other Standard Model processes, including e + e~ — > e + e~ 
and e + e~ — > 77, are not included. Once the primary reaction has been generated, QCD 
radiation and hadronization are done as for hadronic processes. 

e + e~ annihilation to SUSY particles is included as well with complete lowest order dia- 
grams, and cascade decays. The processes include 

e + e~ — > qq 
e + e~ -> E 



e + e~ -> WiWj 
e + e~ — > ZiZj 

e+e" - H L + Z ,H H + Z ,H A + H L ,H A + H° H ,H + + H- 



Note that SUSY Higgs production via WW and ZZ fusion, which can dominate Higgs 
production processes at \fs > 500 GeV, are not included. Also, spin correlations and decay 
matrix elements are not included, as well as initial state photon radiation. 
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e + e cross sections with polarized beams are now included for both Standard Model and 
SUSY processes. The keyword EPOL is used to set Pi(e') and P L (e + ), where 

P L (e) = (n L - n R )/ (n L + n R ) 

so that — 1 < P L < +1. Thus, setting EPOL to —.9, will yield a 95% right polarized electron 
beam scattering on an unpolarized positron beam. 

2.1.10 Technicolor 

Production of a technirho of arbitrary mass and width decaying into W ± Z° or W + W~ pairs. 
The cross section is based on an elastic resonance in the WW cross section with the effective 
W approximation plus a W mixing term taken from EHLQ. Additional technicolor processes 
may be added in the future. 

2.2 QCD Radiative Corrections 

After the primary hard scattering is generated, QCD radiative corrections are added to allow 
the possibility of many jets. This is essential to get the correct event structure, especially at 
high energy. 

Consider the emission of one extra gluon from an initial or a final quark line, 

q(p) -> q(pi)+g(p 2 ) 

From QCD perturbation theory, for small p 2 the cross section is given by the lowest order 
cross section multiplied by a factor 

a = a a s {p 2 )/(27ip 2 )P{z) 

where z = pi/p and P{z) is an Altarelli-Parisi function. The same form holds for the other 
allowed branchings, 

g{p) -> g(pi)+g(p 2 ) 
g(p) -> q(p 1 ) + q(p 2 ) 

These factors represent the collinear singularities of perturbation theory, and they produce 
the leading log QCD scaling violations for the structure functions and the jet fragmentation 
functions. They also determine the shape of a QCD jet, since the jet M 2 is of order a s p 2 
and hence small. 

The branching approximation consists of keeping just these factors which dominate in the 
collinear limit but using exact, non-collinear kinematics. Thus higher order QCD is reduced 
to a classical cascade process, which is easy to implement in a Monte Carlo program. To 
avoid infrared and collinear singularities, each parton in the cascade is required to have a 
mass (spacelike or timelike) greater than some cutoff t c . The assumption is that all physics 
at lower scales is incorporated in the nonperturbative model for hadronization. In ISAJET 
the cutoff is taken to be a rather large value, (6GeV) 2 , because independent fragmentation 
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is used for the jet fragmentation; a low cutoff would give too many hadrons from overlapping 
partons. It turns out that the branching approximation not only incorporates the correct 
scaling violations and jet structure but also reproduces the exact three-jet cross section 
within factors of order 2 over all of phase space. 

This approximation was introduced for final state radiation by Fox and Wolfram. The 
QCD cascade is determined by the probability for going from mass t to mass t 1 emitting 
no resolvable radiation. For a resolution cutoff z c < z < 1 — z c , this is given by a simple 
expression, 

P(t ,t 1 ) = (a s (t )/a s (t 1 )) 2 ^ )/b0 

where 

l{z c ) = C~ ZC dz P(z), b = (33 - 2n/)/(127r) 

J z c 

Clearly if P(t ,ti) is the integral probability, then dP/dti is the probability for the first 
radiation to occur at t\. It is straightforward to generate this distribution and then iteratively 
to correct it to get a cutoff at fixed t c rather than at fixed z c . 

For the initial state it is necessary to take account of the spacelike kinematics and of the 
structure functions. Sjostrand has shown how to do this by starting at the hard scattering 
and evolving backwards, forcing the ordering of the spacelike masses t. The probability 
that a given step does not radiate can be derived from the Altarelli-Parisi equations for the 
structure functions. It has a form somewhat similar to P(t , t\) but involving a ratio of the 
structure functions for the new and old partons. It is possible to find a bound for this ratio 
in each case and so to generate a new t and z as for the final state. Then branchings for 
which the ratio is small are rejected in the usual Monte Carlo fashion. This ratio suppresses 
the radiation of very energetic partons. It also forces the branching g — > t + i for a t quark 
if the t structure function vanishes at small momentum transfer. 

At low energies, the branching of an initial heavy quark into a gluon sometimes fails; 
these events are discarded and a warning is printed. 

Since t c is quite large, the radiation of soft gluons is cut off. To compensate for this, 
equal and opposite transverse boosts are made to the jet system and to the beam jets after 
fragmentation with a mean value 

( Pt 2 ) = (.lGeV)^ 

The dependence on Q 2 is the same as the cutoff used for DRELLYAN and the coefficient is 
adjusted to fit the p t distribution for the W. 

Radiation of gluons from gluinos and scalar quarks is also included in the same approxi- 
mation, but the production of gluino or scalar quark pairs from gluons is ignored. Very little 
radiation is expected for heavy particles produced near threshold. 

Radiation of photons, Ws, and Z's from final state quarks is treated in the same approxi- 
mation as QCD radiation except that the coupling constant is fixed. Initial state electroweak 
radiation is not included; it seems rather unimportant. The W +y s, W~'s and Z's are de- 
cayed into the modes allowed by the WPMODE, WMMODE, and ZOMODE commands respectively. 
Warning: The branching ratios implied by these commands are not included in the cross 
section because an arbitrary number of VF's and Z's can in principle be radiated. 
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2.3 Jet Fragmentation: 



Quarks and gluons are fragmented into hadrons using the independent fragmentation ansatz 
of Field and Feynman. For a quark q, a new quark- ant iquark pair q±qi is generated with 

u : d : s = .43 : .43 : .14 

A meson qq\ is formed carrying a fraction z of the momentum, 

E'+p' z = z(E + p g ) 

and having a transverse momentum p t with (p t ) = 0.35 GeV. Baryons are included by 
generating a diquark with probability 0.10 instead of a quark; adjacent diquarks are not 
allowed, so no exotic mesons are formed. For light quarks z is generated with the splitting 
function 

f(z) = 1 - a + a(b + 1)(1 - z)\ a = 0.96, b = 3 
while for heavy quarks the Peterson form 

f(z) = x(l - x) 2 /((l - xf + ex) 2 

is used with e = .80 /m^ for c and e = .bO/m 2 for q = b,t,y,x. These values of e have 
been determined by fitting PEP, PETRA, and LEP data with ISAJET and should not be 
compared with values from other fits. Hadrons with longitudinal momentum less than zero 
are discarded. The procedure is then iterated for the new quark qi until all the momentum 
is used. A gluon is fragmented like a randomly selected u, d, or s quark or antiquark. 

In the fragmentation of gluinos and scalar quarks, supersymmetric hadrons are not dis- 
tinguished from partons. This should not matter except possibly for very light masses. The 
Peterson form for f(x) is used with the same value of epsilon as for heavy quarks, e = 0.5/m 2 . 

Independent fragmentation correctly describes the fast hadrons in a jet, but it fails to 
conserve energy or flavor exactly. Energy conservation is imposed after the event is generated 
by boosting the hadrons to the appropriate rest frame, rescaling all of the three-momenta, 
and recalculating the energies. 

2.4 Beam Jets 

There is now experimental evidence that beam jets are different in minimum bias events and 
in hard scattering events. ISAJET therefore uses similar a algorithm but different parameters 
in the two cases. 

The standard models of particle production are based on pulling pairs of particles out 
of the vacuum by the QCD confining field, leading naturally to only short-range rapidity 
correlations and to essentially Poisson multiplicity fluctuations. The minimum bias data 
exhibit KNO scaling and long-range correlations. A natural explanation of this was given 
by the model of Abramovskii, Kanchelli and Gribov. In their model the basic amplitude 
is a single cut Pomeron with Poisson fluctuations around an average multiplicity (n), but 
unitarity then produces graphs giving K cut Pomerons with multiplicity K(n). 
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A simplified version of the AKG model is used in ISAJET. The number of cut Pomerons is 
chosen with a distribution adjusted to fit the data. For a minimum bias event this distribution 
is 

P(K) = (1 + AK 2 ) exp -1.8K 

while for hard scattering 

P(l) -> 0.1P(1), P(2) -> 0.2P(2), P(3)^0.5P(3) 

For each side of each event an xq for the leading baryon is selected with a distribution varying 
from flat for K = 1 to like that for mesons for large K: 

f(x) = N(K)(1 - x ) c (K), c(K) = l/K + (1 - l/K)b{s) 

The Xi for the cut Pomerons are generated uniformly and then rescaled to 1 — Xq. Each 
cut Pomeron is then hadronized in its own center of mass using a modified independent 
fragmentation model with an energy dependent splitting function to reproduce the rise in 
dN/dy: 

= 1 _ a + a (b(s) + l)\s), b{s) = b + h log(s) 

The energy dependence is put into f(x) rather than P(K) because in the AKG scheme the 
single particle distribution comes only from the single chain. The probabilities for different 
flavors are taken to be 

u : d : s = .46 : .46 : .08 
to reproduce the experimental K/n ratio. 
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3 Sample Jobs 

The simplest ISAJET job reads a user-supplied parameter file and writes a data file and a 
listing file. The following is an example of a parameter file which generates each type of 
event: 

SAMPLE TWOJET JOB 
800. ,100,2,50/ 
TWOJET 
PT 

50,100,50,100/ 
JETTYPE1 
'GL'/ 
JETTYPE2 

'UP' , 'UB' , 'DN' , 'DB' ,'ST> ,'SB'/ 
END 

SAMPLE DRELLYAN JOB 
800. ,100,2,50/ 
DRELLYAN 
QMW 

80,100/ 
WTYPE 
'W+' , 'W-'/ 
END 

SAMPLE MINBIAS JOB 
800. ,100,2,50/ 
MINBIAS 
END 

SAMPLE WPAIR JOB 
800. ,100,2,50/ 
WPAIR 
PT 

50,100,50,100/ 

JETTYPE1 

'W+' , 'W-' , 'Z0V 

JETTYPE2 

'W+' , 'W-' , 'Z0V 

WM0DE1 

'E+' , 'E-' , 'NUSV 

WM0DE2 

'QUARKS'/ 

END 

SAMPLE HIGGS JOB FOR SSC 

40000,100,1,1/ 

HIGGS 
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QMH 

400,1600/ 

HMASS 

800/ 

JETTYPE1 

'Z0'/ 

JETTYPE2 

'Z07 

WM0DE1 

'MU+' , 'MU-V 

WM0DE2 

'E+' , 'E-V 

PT 

50,20000,50,20000/ 
END 

SAMPLE SUSY JOB 

1800,100,1,10/ 

SUPERSYM 

PT 

50,100,50,100/ 
JETTYPE1 

'GLSS' , 'SQUARKSV 
JETTYPE2 

'GLSS' , ' SQUARKS ' / 
GAUGINO 
60,1,40,40/ 
S QUARK 

80.3,80.3,80.5,81.6,85,110/ 
FORCE 

29,30,1,-1/ 

FORCE 

21,29,1/ 

FORCE 

22,29,2/ 

FORCE 

23,29,3/ 

FORCE 

24,29,4/ 

FORCE 

25,29,5/ 

FORCE 

26,29,6/ 

END 

SAMPLE MSSM JOB FOR TEVATRON 



1800. ,100,1,1/ 

SUPERSYM 

BEAMS 

>P' , 'APV 

MSSMA 

200,-200,500,2/ 
MSSMB 

200,200,200,200,200/ 
MSSMC 

200,200,200,200,200,0,0,0/ 

JETTYPE1 

'GLSS'/ 

JETTYPE2 

'SQUARKS'/ 

PT 

100,300,100,300/ 
END 

SAMPLE MSSM SUGRA JOB FOR LHC 

14000,100,1,10/ 

SUPERSYM 

PT 

50,500,50,500/ 
SUGRA 

247,302,-617.5,10,-1/ 

TMASS 

175/ 

END 

SAMPLE SUGRA HIGGS JOB USING DEFAULT QMH RANGE 

14000,100,20,50/ 

HIGGS 

SUGRA 

200, 200,0, 2, +1/ 

HTYPE 

'HAO'/ 

JETTYPE1 

> GAUGINOS ' , ' SLEPTONS ' / 
JETTYPE2 

> GAUGINOS ' , ' SLEPTONS ' / 
END 

SAMPLE E+E- TO SUGRA JOB 

350. ,100,1,1/ 

E+E- 

SUGRA 

100,100,0,2,-1/ 
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TMASS 

170,-1,1/ 

JETTYPE1 

'ALL'/ 

JETTYPE2 

'ALL'/ 

END 

SAMPLE WH JOB 

2000,100,0,0/ 

WHIGGS 

BEAMS 

' P ' , ' AP ' / 

HMASS 

100./ 

JETTYPE1 

'W+' , 'W-' , 'HTGGS'/ 
JETTYPE2 

'W+' , 'W-' , 'HIGGS'/ 

WM0DE1 

'ALL'/ 

WM0DE2 

'ALL'/ 

PT 

10,300,10,300/ 

END 

STOP 

See Section ^| of this manual for a complete list of the possible commands in a parameter 
file. Note that all input to ISAJET must be in UPPER case only. 

Subroutine RDTAPE is supplied to read events from an ISAJET data file, which is a 
machine-dependent binary file. It restores the event data to the FORTRAN common blocks 
described in Section ^. The skeleton of an analysis job using HBOOK and PAW from the 
CERN Program Library is provided in patch ISAPLT but is not otherwise supported. A 
Zebra output format based on code from the DO Collaboration is also provided in patch 
ISAZEB; see the separate documentation in patch ISZTEXT. 

3.1 DEC VMS 

On a VAX or ALPHA running VMS, ISAJET can be compiled by executing the .COM file 
contained in P=ISAUTIL,D=MAKEVAX. Extract this deck as ISAMAKE.COM and type 

OISAMAKE 

This will run YPATCHY with the pilot patches described in Section f| and the VAX flag to 
extract the source code, decay table, and documentation. The source code is compiled and 
made into a library, which is linked with the following main program, 
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PROGRAM ISARUN 
C MAIN PROGRAM FOR ISAJET 

OPEN (UNIT=1 , STATUS= ' OLD ' , FORM= ' FORMATTED ' READONLY) 

OPEN (UNIT=2 , STATUS= ' NEW , FORM= ' UNFORMATTED ' ) 

OPEN (UNIT=3 , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

OPEN (UNIT=4 , STATUS= ' NEW , FORM= ' FORMATTED ' ) 

CALL ISAJET (-1,2, 3, 4) 

STOP 

END 

to produce ISAJET.EXE. Two other executables, ISASUSY.EXE and ISASUGRA.EXE, 
will also be produced to calculate SUSY masses and decay modes without generating events. 
Temporary files can be removed by typing 

OISAMAKE CLEAN 

Create an input file JOBNAME.PAR following the examples above or the instructions in 
Section ^ and run ISAJET with the command 

OISAJET JOBNAME 

using the ISAJET.COM file contained P=ISAUTIL,D=RUNVAX. This will create a binary 
output file JOBNAME . DAT and a listing file JOBNAME. LIS. Analyze the output data using the 
commands described in Section [8|. 

There is also an simple interactive interface to ISAJET which will prompt the user for 
commands, write a parameter file, and optionally execute it. 

3.2 IBM VM/CMS 

On an IBM mainframe running VM/CMS, run YPATCHY with the pilot patches described 
in Section [| and the IBM flag to extract the source code, decay table, and documentation. 
Compile the source code and link it with the main program 

PROGRAM ISARUN 
C MAIN PROGRAM FOR ISAJET 

OPEN (UNIT=1 , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

OPEN (UNIT=2 , STATUS= ' NEW , FORM= ' UNFORMATTED ' ) 

OPEN (UNIT=3 , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

OPEN (UNIT=4 , STATUS= ' NEW ' , FORM= ' FORMATTED ' ) 

CALL ISAJET (-1,2, 3, 4) 

STOP 

END 

to make ISAJET MODULE. 

Create a file called JOBNAME INPUT containing ISAJET input commands following the ex- 
amples above or the instructions in Section Then run ISAJET using ISAJET EXEC, which 
is contained in P=ISAUTIL,D=RUNIBM. The events will be produced on JOBNAME DATA A 
and the listing on JOBNAME OUTPUT A. 
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3.3 Unix 



The Makefile contained in P=ISAUTIL,D=MAKEUNIX has been tested on DEC Ultrix, 
Hewlett Packard HP-UX, IBM RS/6000 AIX, Linux, Silicon Graphics IRIX, Sun SunOS, 
and Sun Solaris. It should work with minor modifications on almost any Unix system with 
/bin/csh, ypatchy or nypatchy, and a reasonable Fortran 77 compiler. Extract the Makefile 
and edit it, changing the installation parameters to reflect your system. Note in particular 
that CERNlib is usually compiled with underscores postpended to all external names; you 
must choose the appropriate compiler option if you intend to link with it. Then type 

make 

This should produce an executable isajet.x for the event generator, which links the code 
with the following main program: 

PROGRAM RUNJET 

CHARACTER*60 FNAME 

READ 1000, FNAME 
1000 FORMAT (A) 

PRINT 1020, FNAME 
1020 FORMAT (IX, 'Data file = ' ,A) 

OPEN (2 , FILE=FNAME , STATUS= 'NEW , F0RM= ' UNFORMATTED ' ) 

READ 1000, FNAME 

PRINT 1030, FNAME 
1030 FORMAT (IX, 'Parameter file = ' ,A) 

OPEN (3 , FILE=FNAME , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 

READ 1000, FNAME 

PRINT 1040, FNAME 
1040 F0RMAT(1X, 'Listing file = ' ,A) 

OPEN (4 , FILE=FNAME , STATUS= ' NEW , F0RM= ' FORMATTED ' ) 

READ 1000, FNAME 

OPEN ( 1 , FILE=FNAME , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 

CALL ISAJET(-1,2,3,4) 

STOP 

END 

Two other executables, isasusy .x and isasugra.x, will also be produced to calculate SUSY 
masses and decay modes without generating events. Type 

make clean 

to delete the temporary files. 

Most Unix systems do not allow two jobs to read the same decay table file at the same 
time. The shell script in P=ISAUTIL,D=RUNUNIX copies the decay table to a temporary 
file to avoid this problem. Extract this file as isajet. Create an input file jobname.par 
following the examples above or the instructions in Section |G] and run ISAJET with the 
command 
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isajet jobname 



This will create a binary output file jobname.dat and a listing file jobname. lis. Analyze 
the output data using the commands described in Section || 

This section only describes running ISAJET as a standalone program and generating 
output in machine-dependent binary form. The user may elect to analyze events as they are 
generated; this is discussed in Section |5| of this manual. 
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4 Patchy and PAM Organization 



Patchy is a code management system developed at CERN and used to maintain the CERN 
Library. It is used to provide versions of ISAJET for a wide variety of computers. Instructions 
for using PATCHY are available from http://wwwinfo.cern.ch/asdoc/Welcome.html. 

A master source file in Patchy is called a "PAM." The ISAJET PAM contains all the 
source code and documentation plus Patchy commands to include common blocks and to 
select the desired version. It is divided into the following patches: 

ISACDE: contains all common blocks, etc. These are divided into decks based on their 
usage. 

ISADATA: contains block data ALDATA. This must always be loaded when using ISAJET. 
ISAJET: contains the code for generating events. Each subroutine is in a separate deck 
with the same name. 

ISASSRUN: contains the main program for IS ASUS Y, which prompts for input parameters 
and prints out all the decay modes. It is selected by * ISASUSY. 

ISASUSY: contains code to calculate all the decay widths and branching fractions in the 
minimal supersymmetric model. 

I SAT APE: contains the code for reading and writing tapes, again with each subroutine on 
a separate deck. 

ISARUN: contains a main program and a simple interactive interface. It is selected by 
IF=INTERACT. 

ISAZEB: contains Zebra format output routines, an alternative to the ISATAPE routines. 
ISZRUN: contains the analog of ISAPLT for the Zebra format. 

ISAPLT: contains a simple calorimeter simulation and the skeleton of a histogramming 
job using HBOOK. 

ISATEXT: contains the instructions for using ISAJET, i.e. the text of this document. It 
also includes the documentation for ISASUSY. 

ISZTEXT: contains the instructions for the Zebra output routines and a description of the 
Zebra banks. 

ISADECAY: contains the input decay table. 

The code is actually maintained using RCS on a Silicon Graphics computer at BNL. 
Patchy is used primarily to handle common blocks and machine dependent code. 

The input to YPATCHY must contain both +USE cards, which define the wanted program 
version, and +EXE cards, which determine which patches or decks are written to the ASM 
file. To facilitate this selection, the ISAJET PAM contains the following pilot patches: 

* ISADECAY: USE selects ISADECAY and all corrections to it. 

*ISAJET: USE selects ISACDE, ISADATA, ISAJET, ISATAPE, ISARUN and all correc- 
tions to them. Note that ISARUN is not actually selected without +USE, INTERACT. 

*ISAPLT: USE selects ISACDE, ISAPLT, and all corrections to them. 

*ISASUSY: USE selects CDESUSY, ISASUSY, and ISASSRUN to create a program to 
calculate all the MSSM decay modes. 

*ISATEXT: USE selects ISACDE, ISATEXT, and all corrections to them. 

*ISAZEB: USE selects ISAJET with a Zebra output format. 

*ISZRUN: USE selects the Zebra analysis package. 
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Patches are provided to select the machine dependent features for specific computers or 
operating systems: 

ANSI: ANSI standard Fortran (no time or date functions) 

APOLLO: APOLLO - only tested by CERN 

CDC: CDC 7600 and 60-bit CYBER (obsolete) 

CRAY: CRAY with UNICOS 

DECS: DEC Station with Ultrix 

ETA: ETA 10 running Unix System V (obsolete) 

HPUX: HP/9000 7xx running Unix System V 

IBM: IBM 370 and 30xx running VM/CMS 

IBMRT: IBM RS/6000 running AIX 3.x or 4.x 

LINUX: PC running Linux with f2c/gcc or g77 compiler 

SGI: Silicon Graphics running IRIX 

SUN: Sun Sparcstation running SUNOS or Solaris 

VAX: DEC VAX or Alpha running VMS 
These patches in turn select a variety of patches and IF flags, allowing one to select more 
specific features, as indicated below. (Replace & by + everywhere.) 



WATCH, ANSI. 

&USE, DOUBLE. 

&USE,STDI0. 

&USE,M0VEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE , RANFCALL . 

&USE , NOCERN , IF=-CERN . 

&E0D 



GENERIC ANSI FORTRAN. 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 



&PATCH, APOLLO . 

&DECK , BLANKDEK . 

&USE, DOUBLE. 

&USE,STDI0. 

&USE,M0VEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&USE,IMPN0NE. 

&E0D. 



DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



&PATCH,CDC. 

&USE, SINGLE. 

&USE,LEVEL2. 

&USE,CDCPACK. 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&E0D 



CDC 7600 OR CYBER 175. 
SINGLE PRECISION. 
LEVEL 2 STORAGE. 

PACK 2 WORDS PER WORD FOR INPUT/OUTPUT. 
STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 



24 



&PATCH , CRAY . 

&USE, SINGLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH,DECS. 

&USE,SUN. 

&EOD 

&PATCH , ETA . 

&USE, SINGLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFCALL . 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH,HPUX. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 

&PATCH,IBM. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH , IBMRT . 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 



CRAY XMP OR 2. 
SINGLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
NO CERN LIBRARY. 



DEC STATION (ULTRIX) 



ETA- 10. 

SINGLE PRECISION. 
STANDARD FORTRAN 77 
FORTRAN REPLACEMENT 
STANDARD RANSET AND 
NO CERN LIBRARY. 



TAPE INPUT/ OUTPUT. 
FOR MOVLEV. 
RANGET CALLS. 



HP/9000 7XX RUNNING UNIX. 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



IBM 370 OR 30XX. 
DOUBLE PRECISION. 
STANDARD FORTRAN 77 
FORTRAN REPLACEMENT 
FORTRAN RANF. 
STANDARD RANSET AND 
NO CERN LIBRARY. 



TAPE INPUT/OUTPUT. 
FOR MOVLEV. 

RANGET CALLS. 



IBM RS/6000 WITH AIX 3.1 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 
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&USE , RANFCALL . STANDARD RANSET AND RANGET CALLS. 

&USE, NOCERN, IF=-CERN. NO CERN LIBRARY. 

&USE,IMPNONE. IMPLICIT NONE 
&EOD 



&PATCH, LINUX. 
&USE, DOUBLE. 

&USE,STDIO. STANDARD FORTRAN 77 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 



IBM PC WITH LINUX l.X 
DOUBLE PRECISION. 
TAPE INPUT/ OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



&PATCH.SGI. 

SILICON GRAPHICS 4D/XX. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&EOD 



DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 



&PATCH,SUN. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&EOD 



SUN (SPARC) 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 



&PATCH,VAX. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 



DEC VAX 11/780 OR 8600. 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



An empty patch INTERACT selects a main program and an interactive interface which 
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will prompt the user for parameters and do some error checking. A patch CERN allows 
ISAJET to take the random number generator RANF and several other routines from the 
CERN Library; to use this include the Patchy command 

&USE,CERN. 

Similarly, a patch PDFLIB enables the interface to the PDFLIB parton distribution compi- 
lation by H. Plothow-Besch: 

&USE, PDFLIB 

The only internal links with PDFLIB are calls to the routines PDFSET, PFTOPDG, and 
DXPDF, and the common blocks W50510 and W50517, 

C Copy of PDFLIB common block 

C0MM0N/W50510/IFLPRT 

INTEGER IFLPRT 

SAVE /W50510/ 
C Copy of PDFLIB common block 

C0MM0N/W50517/N6 

INTEGER N6 

SAVE /W50517/ 

which are used to specify the level of output messages and the logical unit number for them. 

In general it should be sufficient to run YPATCHY with the following cradle (replace & 
with + everywhere): 

&USE, (*ISAJET, *ISATEXT, *ISADECAY, *ISAPLT) . CHOOSE ONE. 

&USE, ANSI, DECS, HPUX, IBM, IBMRT, SGI, SUN, CHOOSE ONE. 

& [USE, INTERACT] . FOR INTERACTIVE MODE. 

& [USE, CERN.] FOR CERN LIBRARY. 

&[USE,HB00K3.] HBOOK 3 FOR ISAPLT. 

&EXE. 

&PAM. 

&QUIT. 

The input to YPATCHY can also contain changes by the user. It is suggested that these 
not be made permanent parts of the PAM to avoid possible conflicts with later corrections. 
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5 Main Program 



A main program is not supplied with ISAJET. To generate events and write them to disk, 
the user should provide a main program which opens the files and then calls subroutine 
ISAJET. In the following sample, i,j,m,n are arbitrary unit numbers. 
Main program for VMS: 

PROGRAM RUNJET 

C 

C MAIN PROGRAM FOR ISAJET ON BNL VAX CLUSTER. 

C 

0PEN(UNIT=i,FILE='$2$DUA14: [ISAJET . ISALIBRARY] DECAY . DAT ' , 
$STATUS= ' OLD ' , F0RM= ' FORMATTED ' , READONLY) 
OPEN (UNIT= j , FILE= ' my j ob . dat ' , STATUS= ' NEW ' , F0RM= ' UNFORMATTED ' ) 
OPEN (UNIT=m , FILE= ' my j ob . par ' , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 
OPEN (UNIT=n , FILE= ' my j ob . 1 i s ' , STATUS= ' NEW ' , F0RM= ' FORMATTED ' ) 

C 

CALL ISAJET(+-i,+-j ,m,n) 

C 

STOP 
END 

Main program for IBM (VM/CMS) 
PROGRAM RUNJET 

C 

C MAIN PROGRAM FOR ISAJET ON IBM ASSUMING FILES HAVE BEEN 

C OPENED WITH FILEDEF. 

C 

CALL ISAJET (+-i ,+- j ,m,n) 

C 

STOP 
END 

Main program for Unix: 
PROGRAM RUNJET 

C 

C Main program for ISAJET on Unix 

C 

CHARACTER*60 FNAME 

C 

C Open user files 

READ 1000, FNAME 
1000 FORMAT (A) 

PRINT 1020, FNAME 
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1020 FORMAT (IX, 'Data file = ' ,A) 

OPEN (2 , FILE=FNAME , STATUS= ' NEW ' , FORM= ' UNFORMATTED ' ) 

READ 1000, FNAME 

PRINT 1030, FNAME 
1030 FORMAT (IX, 'Parameter file = ' ,A) 

OPEN (3 , FILE=FNAME , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

READ 1000, FNAME 

PRINT 1040, FNAME 
1040 FORMAT (IX, 'Listing file = ' ,A) 

OPEN (4 , FILE=FNAME , STATUS= 'NEW, FORM= ' FORMATTED ' ) 
C Open decay table 

READ 1000, FNAME 

OPEN ( 1 , FILE=FNAME , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

C 

C Run ISA JET 

CALL ISAJET(-1,2,3,4) 

C 

STOP 
END 

The arguments of ISAJET are tape numbers for files, all of which should be opened by 
the main program. 

TAPEi: Decay table (formatted). A positive sign prints the decay table on the output 
listing. A negative sign suppress printing of the decay table. 

TAPEj: Output file for events (unformatted). A positive sign writes out both resonances 
and stable particles. A negative sign writes out only stable particles. 

TAPEm: Commands as defined in Section 6 (formatted). 

TAPEn: Output listing (formatted). 
In the sample jobs in Section 3, TAPEm is the default Fortran input, and TAPEn is the 
default Fortran output. 

5.1 Interactive Interface 

To use the interactive interface, replace the call to ISAJET in the above main program by 

CALL ISASET (+-i, +- j ,m,n) 
CALL ISAJET(+-i,+-j ,m,n) 

ISASET calls DIALOG, which prompts the user for possible commands, does a limited 
amount of error checking, and writes a command file on TAPEm. This command file is 
rewound for execution by ISAJET. A main program is included in patch ISARUN to open 
the necessary files and to call ISASET and ISAJET. 
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5.2 User Control of Event Loop 



If the user wishes to integrate ISAJET with another program and have control over the event 
generation, he can call the driving subroutines himself. The driving subroutines are: 

ISAINI(+-i,+-j ,m,n): initialize ISAJET. The arguments are the same as for subroutine 
ISAJET. 

ISABEG(IFL): begin a run. IFL is a return flag: IFL=0 for a good set of commands; 
IFL=1001 for a STOP; any other value means an error. 

ISAEVT(I , OK, DONE) generate event I. Logical flag OK signifies a good event (almost 
always .TRUE.); logical flag DONE signifies the end of a run. 

ISAEND: end a run. 

There are also subroutines provided to write standard ISAJET records, or Zebra records if 
the Zebra option is selected: 

ISAWBG to write a begin-of-run record, should be called immediately after ISABEG 
ISAWEV to write an event record, should be called immediately after ISAEVT 
ISAWND to write an end-of-run record, should be called immediately after ISAEND 
The control of the event loop is somewhat complicated to accomodate multiple evolu- 
tion and fragmentation as described in Section 11. Note in particular that after calling 
ISAEVT one should process or write out the event only if OK=.TRUE. The check on the 
DONE flag is essential if one is doing multiple evolution and fragmentation. The following 
example indicates how events might be generated, analyzed, and discarded (replace & by + 
everywhere) : 

PROGRAM SAMPLE 

C 

&SELF,IF=IMPNONE 

IMPLICIT NONE 

&SELF 

&CDE,ITAPES 

&CDE,IDRUN 

&CDE,PRIMAR 

&CDE,ISL00P 

C 

INTEGER JTDKY , JTEVT , JTCOM , JTLIS , IFL , IL00P 
LOGICAL OK, DONE 
SAVE IL00P 

c 

C> Open files as above 

C> Call user initialization 

c 

C 

C Initialize ISAJET 

C 

CALL ISAINI(-i,0,m,n) 
1 IFL=0 
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CALL ISABEG(IFL) 
IF(IFL.NE.O) STOP 

C 

C Event loop 

C 

IL00P=0 
101 CONTINUE 

IL00P=IL00P+1 
C Generate one event - discard if .NOT. OK 

CALL ISAEVT(IL00P,0K,D0NE) 
IF (OK) THEN 

c 

C> Call user analysis for event 

c 

END IF 

IF (.NOT. DONE) GO TO 101 

C 

C Calculate cross section and luminosity 

C 

CALL ISAEND 

c 

C> Call user summary 

c 

GO TO 1 
END 

5.3 Multiple Event Streams 

It may be desirable to generate several different kinds of events simultaneously to study 
pileup effects. While normally one would want to superimpose minimum bias or low-pt jet 
events on a signal of interest, other combinations might also be interesting. It would be very 
inefficient to reinitialize IS A JET for each event. Therefore, a pair of subroutines is provided 
to save and restore the context, i.e. all of the initialization information, in an array. The 
syntax is 

CALL CTX0UT(NC,VC,MC) 
CALL CTXIN(NC,VC,MC) 

where VC is a real array of dimension MC and NC is the number of words used, about 20000 
in the standard case. If NC exceeds MC, a warning is printed and the job is terminated. The 
use of these routines is illustrated in the following example, which opens the files with names 
read from the standard input and then superimposes on each event of the signal sample 
three events of a pileup sample. It is assumed that a large number of events is specified in 
the parameter file for the pileup sample so that it does not terminate. 

PROGRAM SAMPLE 



31 



c 

C Example of generating two kinds of events. 

C 

CHARACTER*60 FNAME 
REAL VC1 (20000) ,VC2 (20000) 
LOGICAL 0K1,D0NE1,0K2,D0NE2 
INTEGER NCI , NC2 , IFL , IL00P ,12, IL00P2 

C 

C Open decay table 

READ 1000, FNAME 
1000 FORMAT (A) 

OPEN ( 1 , FILE=FNAME , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 
C Open user files 

READ 1000, FNAME 

OPEN (3 , FILE=FNAME , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 
READ 1000, FNAME 

OPEN (4 , FILE=FNAME , STATUS= 'NEW , FORM= ' FORMATTED ' ) 
READ 1000, FNAME 

OPEN ( 13 , FILE=FNAME , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 
READ 1000, FNAME 

OPEN ( 14 , FILE=FNAME , STATUS= 'NEW , FORM= ' FORMATTED ' ) 

C 

C Initialize ISAJET 

CALL ISAINK-1,0,3,4) 
CALL CTXOUT (NCI, VC1, 20000) 
CALL ISAINK-1,0,13,14) 
IFL=0 

CALL ISABEG(IFL) 

IF(IFL.NE.O) ST0P1 

CALL CTXOUT (NC2,VC2, 20000) 

IL00P2=0 

CALL user_initialization_routine 

C 

1 IFL=0 

CALL CTXIN(NC1,VC1, 20000) 

CALL ISABEG(IFL) 

CALL CTXOUT (NCl.VC 1,20000) 

IF(IFL.NE.O) GO TO 999 

ILOOP=0 

C 

C Main event 

C 

101 CONTINUE 

IL00P=IL00P+1 
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CALL CTXIN (NCI ,VC1, 20000) 
CALL ISAEVT(IL00P,0K1,D0NE1) 
CALL CTXOUT (NCI, VC1, 20000) 
IF(.NOT.OKl) GO TO 101 
CALL user_analysis_routine 

C 

C Pileup 
C 

CALL CTXIN (NC2,VC2, 20000) 
12=0 

201 CONTINUE 

IL00P2=IL00P2+1 

CALL ISAEVT(IL00P2,0K2,D0NE2) 

IFC0K2) 12=12+1 

IF(D0NE2) ST0P2 

CALL user_analysis_routine 
IFCI2.LT. 3) GO TO 201 
CALL CTXOUT (NC2,VC2, 20000) 

C 

IF(.NOT.DONEl) GO TO 101 

C 

C Calculate cross section and luminosity 

C 

CALL CTXIN (NCI ,VC1 ,20000) 
CALL ISAEND 
GO TO 1 

C 

999 CALL CTXIN (NC2,VC2, 20000) 
CALL ISAEND 

CALL user_termination_routine 

STOP 

END 

It is possible to superimpose arbitrary combinations of events, including events of the same 
reaction type with different parameters. In general the number of events would be selected 
randomly based on the cross sections and the luminosity. 

At this time CTXOUT and CTXIN cannot be used with the Zebra output routines. 
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6 Input 



ISAJET is controlled by commands read from the specified input file by subroutine READIN. 
(In the interactive version, this file is first created by subroutine DIALOG.) Syntax errors 
will generate a message and stop execution. Based on these commands, subroutine LOGIC 
will setup limits for all variables and check for inconsistencies. Several runs with different 
parameters can be combined into one job. The required input format is: 

Title 

Ecm , Nevent , Npr int , N j ump/ 
Reaction 

(Optional parameters) 
END 

(Optional additional runs) 
STOP 

with all lines starting in column 1 and typed in upper case. These lines are explained below. 

Title line: Up to 80 characters long. If the first four letters are STOP, control is returned 
to main program. If the first four letters are SAME, the parameters from previous run are 
used excepting those which are explicitly changed. 

Ecm line: This line must always be given even if the title is SAME. It must give the 
center of mass energy (Ecm) and the number of events (Nevent) to be generated. One may 
also specify the number of events to be printed (Nprint) and the increment (Njump) for 
printing. The first event is always printed if Nprint > 0. For example: 

800. ,1000,10,100/ 

generates 1000 events at E cm = 800 GeV and prints 10 events. The events printed are: 
1,100,200,. . . . Note that an event typically takes several pages of output. This line is read 
with a list directed format (READ*). 

After Nprint events have been printed, a single line containing the run number, the event 
number, and the random number seed is printed every Njump events (if Njump is nonzero). 
This seed can be used to start a new job with the given event if in the new run NSIGMA is 
set equal to zero: 

SEED 
value/ 
NSIGMA 
0/ 

In general the same events will only be generated on the same type of computer. 

Reaction line: This line must be given unless title is SAME, when it must be omitted. 
It selects the type of events to be generated. The present version can generate TWOJET, 
E+E-, DRELLYAN, MINBIAS, WPAIR, SUPERSYM, HIGGS, PHOTON, TCOLOR, or 
WHIGGS events. This line is read with an A8 format. 

Optional parameters: Each optional parameter requires two lines. The first line is a 
keyword specifying the parameter and the second line gives the values for the parameter. 
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The parameters can be given in any order. Numerical values are read with a list directed 
format (READ*), jet and particle types are read with a character format and must be 
enclosed in quotes, and logical flags with an LI format. All momenta are in GeV and all 
angles are in radians. 

The parameters can be classified in several groups: 



Jet Limits: 


W/H Limits: 


Decays: 


Constants: 


Other: 


JETTYPE1 


HTYPE 


FORCE 

1 V ./ 1 V V. A i 


CUTJET 

v — V... 1 \J A i 1 


BEAMS 


JETTYPE2 

l_} A i _l_ JL J_ 1 A j 


PHIW 


FORCE1 

_L V ./ A. v v. A i 1 


CUTOFF 


EPOL 


JETTYPE3 

l_} A i .1 J- 1 1 A iKJ 


QMH 


NODECAY 

1 1 V. ./ A s A i * — ± 1 1 


FRAGMENT 

1 1 l.. 1 V — \ ill J i 1 1 _l_ 


NPOMERON 

111 V ./ 1 V 1 J i 1 \ V .X 1 1 


P 


QMW 


NOETA 


GAUGINO 

V_ J ± 1 \_J V — \ 11 1 V. ./ 


NSIGMA 

1 1 K-J 1 * — ' IVli 1 


PHI 


OTW 


NOEVOLVE 


GMSB 


NTRIES 


PT 


THW 


NOFRGMNT 

1 1 V ./ 1 A. V * — \ Ull 1 _J_ 


GMSB2 

V — \ ill K-J A I ^ 


PDFLIB 


TH 


WTYPE 


NOGRAV 

1 1 V ./ V — \ 1 \j± 1 V 


HMASS 

lll\li 1UU 


SEED 


x 


XW 


NOPIO 


HMASSES 

lllVli IkJkJ J — / 


STRUC 


Y 


YW 




LAMBDA 


WFUDGE 

v v i vy j — ' \_» j — ; 


WMODE1 






MGVTNO 


WMMODE 


WMODE2 

TV 1V1 V ./ A s A i ^ 






MSSMA 

lllUkJIVli 1 


WPMODE 

V V 1 1 V 1 V ./ A s A i 








MSSMB 

1V1 K-J K^J 1 V 1 1 J 


ZOMODE 

i—l \J 1 V 1 V ./ A s A i 








MSSMC 










MSSMD 

1V1U kJ 1 VI J / 










MSSME 

1 V 1 K-J K-J 1 V 1 1 i 










NUSUG1 










NUSUG2 

1 'I KJ K-> KJ vJ 










NUSUG3 










NUSUG4 










NUSUG5 










SIGQT 










SIN2W 










SLEPTON 










SQUARK 










SUGRA 










TCMASS 










TMASS 










WMASS 










XGEN 





It may be helpful to know that the TWOJET, WPAIR, PHOTON, SUPERSYM, and 
WHIGGS processes use the same controlling routines and so share many of the same vari- 
ables. In particular, PT limits should normally be set for these processes, and JETTYPE1 
and JETTYPE2 are used to select the reactions. Similarly, the DRELLYAN, HIGGS, and 
TCOLOR processes use the same controlling routines since they all generate s-channel reso- 
nances. The mass limits for these processes are set by QMW. Normally the QMW limits will 
surround the W ± , Z°, or Higgs mass, but this is not required. (QMH acts like QMW for the 
Higgs process.) For historical reasons, JETTYPE1 and JETTYPE2 are used to select the W 
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decay modes in DRELLYAN, while WMODE1 and WMODE2 select the W decay modes for 
WPAIR, HIGGS, and WHIGGS. Also, QTW can be used to generate DRELLYAN events 
with non-zero transverse momentum, whereas HIGGS automatically fixes QTW to be zero. 
(Of course, non-zero transverse momentum will be generated by gluon radiation.) 
A complete list of keywords and their default values follows. 
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Keyword 
Values 


Default values 


Explanation 


BEAMS 
typei,type 2 


'P','P' 


Initial beams. Allowed are 
'P','AP','N','AN'. 


CUT JET 

He 


6. 


Cutoff mass for QCD jet 
evolution. 


CUTOFF 


.200,1.0 


Cutoff qt 2 = ij 2 Q u for 
DRELLYAN events. 


EPOL 

P L (e-),P L (e+) 


0,0 


Polarization of e~ (e + ) beam, 
P L (e) = (n L -n R )/(n L -n R ), 
so that -1 < P L < 1 


FORCE 

Mi, —,hl 


None 


Force decay of particles, 
±i -> ±(il + ... +i5). 
Can call 20 times. 
See note for i = quark. 


FORCE1 

Mi, -,15/ 


None 


Force decay % — > il + ... + i5. 
Can call 40 times. 
See note for % = quark. 


FRAGMENT 

Pudi- ■ ■ 


.4,... 


Fragmentation parameters. 
See also SIGQT, etc. 


GAUGINO 

mi,m 2 ,m3,m 4 


50,0,100,100 


Masses for g, 7, 
W + , and Z° 


GMSB 

A m ,MmM 

tan /3,sgn /i,C gr 


none 


GMSB messenger SUSY breaking, 
mass, number of 5 + 5, VEV 
ratio, sign, gravitino scale 


GMSB2 

ft,5M 2 H JM z Hu ,D Y {M) 
N 5l ,N 52 ,N 5s 


1 n n ri 

1,0,0,0 


non-minimal GMSB parameters 

1 j_ * 1 * 

gaugmo mass multiplier 

Higgs mass shifts, D-term mass 2 

indep. gauge group messengers 


HMASS 

m 





Mass for standard Higgs. 


HMASSES 
mi,. . . ,m 9 


0,...,0 


Higgs meson masses for 
charges 0,0,0,0,0,1,1,2,2. 
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HTYPE 
'HLO'/ or... 

JETTYPE1 
'GL\'UP\... 

JETTYPE2 
'GL\'UP\... 

JETTYPE3 
'GL','UP',... 



none 



'ALL' 



'ALL' 



'ALL' 



One MSSM Higgs type ('HLO', 
'HHO', or 'HAO') 

Select types for jets: 
'ALL'; 'GL'; 'QUARKS'='UP', 
'UB','DN','DB','ST','SB', 
'CH','CB','BT','BB','TP', 
'TB','X','XB','Y','YB'; 
'LEPTONS'='E-','E+','MU-', 
'MU+','TAU-','TAU+'; 'NUS'; 
'GM','W+','W-','Z0' 
See note for SUSY types. 



LAMBDA 
A 



QCD scale 



MGVTNO 

-^gravitino 

MSSMA 
m(g),n, 
m(A), tan (3 

MSSMB 

m(qi) ,m(d r ) ,m(u r ) , 
m(li),m(e r ) 

MSSMC 

m(q 3 ),m(b r ),m(t r ), 
m(lz) ,m(r r ) , 

MSSMD 

m(q 2 ),m(s r ),m(c r ), 
mfc) ,m(mu r ) 

MSSME 
M..M-. 



10 20 GeV 



None 



None 



None 



Gravitino mass - ignored for 
GMSB model 

MSSM parameters - required 
Gluino mass, /i, A mass, 
tan (5 

MSSM 1st generation - required 
Left and right soft squark and 
slepton masses 

MSSM 3rd generation - required 
Soft squark masses, slepton 
masses, and squark and slepton 
mixings 



MSSM 2nd generation - optional 
from MSSMB Left and right soft squark and 
slepton masses 

MSSM gaugino masses - optional 
GUT values Default is to scale from gluino 



NODECAY 
TRUE or FALSE 



FALSE 



Suppress all decays. 



NOETA 
TRUE or 



FALSE FALSE 



Suppress eta decays. 
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NOEVOLVE 

TRUE or FALSE FALSE 
NOGRAV 

TRUE or FALSE FALSE 
NOHADRON 

TRUE or FALSE FALSE 
NONUNU 

TRUE or FALSE FALSE 
NOPIO 

TRUE or FALSE FALSE 

NPOMERON 

ni,n 2 1,20 

NSIGMA 

n 20 
NTRIES 

n 1000 
NUSUG1 

M 1 ,M 2 ,M 3 none 
NUSUG2 

A t ,Ab,A T none 
NUSUG3 

M Hd ,M Hu none 
NUSUG4 

M UL ,M dR ,M UR , none 
NUSUG5 

M tL ,M bR ,M tR , none 
M TL ,M TR 

P 

Pmin(l)r • • ;Pmax(3) l.,0.5E cm 



Suppress QCD evolution and 
hadronization. 

Suppress gravitino decays in 
GMSB model 

Suppress hadronization of 
jets and beam jets. 

Suppress Z° neutrino decays. 

Suppress it decays. 

Allow rii < n < n 2 cut pomerons. 
Controls beam jet mult. 

Generate n unevolved events 
for SIGF calculation. 

Stop if after n tries 
cannot find a good event. 

Optional non-universal SUGRA 
gaugino masses 

Optional non-universal SUGRA 
A terms 

Optional non-universal SUGRA 
Higgs masses 

Optional non-universal SUGRA 
1st /2nd generation masses 

Optional non-universal SUGRA 
3rd generation masses 

Momentum limits for jets. 
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PDFLIB 

'namei',vali,. . . None 
PHI 

0min(l),- • • ,0max(3) 0,27T 

PHIW 

0min,0max 0,27T 

PT or PPERP 

Pt,mm 

(3) .05£' cm ,.2£' cm 

QMH 

Q'miiuQ'max 

QMW 

QmimQmax 

QTW 

Qt, 111111;%, max 

SEED 
real 

SIGQT 
a 

SIN2W 
sin 2 (0 w ) 

SLEPTON 
mi,. . . ,m 6 

SQUARK 
mi,. . . ,m 6 

STRUC 
name 

SUGRA 

m ,mi/ 2 ,Ao, none 
tan /9,sgn /i 



CERN PDFLIB parton distribution 
parameters. See PDFLIB manual. 

Phi limits for jets. 
Phi limits for W. 



Pt limits for jets. 

Default for TWOJET only. 



Minimal supergravity parameters 
scalar M, gaugino M, trilinear 
breaking term, vev ratio, +-1 



Mass limits for Higgs. 
.05£ cm ,.2£ cm Equivalent to QMW. 

Mass limits for W. 

.05E cm ,.2E cm 

q t limits for W. Fix q t = 
.l,.025£' cm for standard Drell-Yan. 

Random number seed (double 
precision if 32 bit). 

Internal k t parameter for 
.35 jet fragmentation. 



.232 



Weinberg angle. See WMASS. 



Masses for u e , e, 
100,. ..,101.8 /j, i> T , f 

Masses for u, d, s, 
100.3,. ..,240. c, 6, t 

Structure functions. CTEQ3L, 
'CTEQ3L' CTEQ2L, EHLQ, OR DO 
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TH or THETA 

0min(l)v • • ,0max(3) 0,7T 

THW 

^min^max 0,7T 

TCMASS 

m,T 1000,100 
TMASS 

m t ,m y ,m x 180.,-1.,-1. 
WFUDGE 

factor 1.85 



Theta limits for jets. Do not 
also set Y. 

Theta limits for W. Do not 
also set YW. 

Technicolor mass and width. 

t, y, and x quark masses. 

Fudge factor for DRELLYAN 
evolution scale. 



WMASS 

M W ,M Z 80.2, 91.19 

WMMODE 

'UP',. . . ,'TAU+' 'ALL' 
WMODE1 

'UP7UB',... 'ALL' 
WMODE2 

'UP','UB',. .. 'ALL' 
WPMODE 

'UP',. . . ,'TAU+' 'ALL' 
WTYPE 

typei.typea 'GM','Z0' 



W and Z masses. See SIN2W. 



Decay modes for W in parton 
cascade. See JETTYPE. 

) 

) Decay modes for WPAIR. 
)Same code for quarks and 
)leptons as JETTYPE. 

) 

Decay modes for W + in parton 
cascade. See JETTYPE. 

Select W type: W+,W-,GM,Z0. 
Do not mix W+,W- and GM,Z0. 



X Feynman x limits for jets. 

■£min(l)v • • )-£max(3) 1,1 



XGEN 

a(l),...,a(8) .96,3,0,.8,.5,... 
XGENSS 

a(l),...,a(7) .5,-5,. .. 



Jet fragmentation, Peterson 
with e = a(n)/m 2 , n = 4-8. 

Fragmentation of GLSS, UPSS, 
etc. with e = a{n)/m * *2 
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XW Feynman x limits for W. 

Y Y limits for each jet. 

ymin(l),- • • ,y m ax(3) from PT Do not also set TH. 

YW Y limits for W. 

2/min,2/max from QTW,QMW Do not set both YW and THW. 

ZOMODE Decay modes for Z° in parton 

'UP',. . . ,'TAU+' 'ALL' cascade. See JETTYPE. 
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For example the lines 

P 

40. ,50. ,10. ,100./ 

would set limits for the momentum of jet 1 between 40 and 50 GeV, and for jet 2 between 
10 and 100 GeV. As another example the lines 

WTYPE 
'W+'/ 

would specify that for DRELLYAN events only W+ events will be generated. If for a 
kinematic variable only the lower limit is specified then that parameter is fixed to the given 
value. Thus the lines 

P 

40. ,,10./ 

will fix the momentum for jet 1 to be 40 GeV and for jet 2 to be 10 GeV. If only the upper 
limit is specified then the default value is used for the lower limit. Jet 1 or jet 2 parameters 
for DRELLYAN events refer to the W decay products and cannot be fixed. If QTW is fixed 
to 0, then standard Drell-Yan events are generated. 

Note that if the limits given cover too large a kinematic range, the program can become 
very inefficient, since it makes a fit to the cross section over the specified range. NTRIES 
has to be increased if narrow limits are set for X, XW or for jet 1 and jet 2 parameters 
in DRELLYAN events. For larger ranges several runs can be combined together using the 
integrated cross section per event SIGF/NEVENT as the weight. This cross section is 
calculated for each run by Monte Carlo integration over the specified kinematic limits and 
is printed at the end of the run. It is corrected for JETTYPEi, WTYPE, and WMODEi 
selections; it cannot be corrected for branching ratios of forced decays or for WPMODE, 
WMMODE, or Z0MODE selections, since these can affect an arbitrary number of particles. 

For e+e" -> W + W~, Z°Z°, use FORCE and FORCE1 instead of WMODEi to select 
the W decay modes. Note that these do not change the calculated cross section. (In the 
E+E- process, the W and Z decays are currently treated as particle decays, whereas in the 
WPAIR and HIGGS processes they are treated as 2 — > 4 parton processes.) 

Normally the user should set PT limits for TWOJET, PHOTON, WPAIR, SUPERSYM, 
and WHIGGS events and QMW and QTW limits for DRELLYAN, HIGGS, and TCOLOR 
events. If these limits are not set, they will be selected as fractions of E cm . This can give 
nonsense. For TWOJET the p t range should usually be less than about a factor of two 
except for b and t jets at low p t to produce uniform statistics. For W + , W~, or Z° events 
or for Higgs events the QMW (QMH) range should usually include the mass. But one can 
select different limits to study, e.g., virtual W production or the effect of a lighter or heavier 
Higgs on WW scattering. If only t decays are selected, then the lower QMW limit must be 
above the t threshold. For standard Drell-Yan events QTW should be fixed to zero, 

QTW 

0/ 
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Transverse momenta will then be generated by initial state gluon radiation. A range of QTW 
can also be given. For SUPERSYM either the masses and decay modes should be specified, 
or the MSSM, SUGRA, or GMSB parameters should be given. For fourth generation quarks 
it is necessary to specify the quark masses. 

SAME cannot be used to combine standard DRELLYAN events (QTW fixed equal to 0) 
and DRELLYAN events with nonzero QTW. 

For HIGGS with W + W~ or Z°Z° decays allowed it is generally necessary to set PT limits 
for the Ws, e.g. 

PT 

50,20000,50,20000/ 

If this is not done, then the default lower limit of 1 GeV is used, and the t-channel exchanges 
will dominate, as they should in the effective W approximation. Depending on the other 
parameters, the program may fail to generate an event in NTRIES tries. 

SUPERSYM (SUSY) by default generates just gluinos and squarks in pairs. There are no 
default masses or decay modes. Masses can be set using GAUGINO, SQUARE, SLEPTON, 
and HMASSES. Decay modes can be specified with FORCE or by modifying the decay 
table. Left and right squarks are distinguished but assumed to be degenerate, except for 
stops. Since version 7.11, types must be selected with JETTYPEi using the supersymmetric 
names, e.g. 

JETTYPEI 

'GLSS' , 'UPSSL' , 'UPSSRV 

Use of the corresponding standard model names, e.g. 

JETTYPEI 
'GL' , 'UP'/ 

and generation of pure photinos, winos, and zinos are no longer supported. 

If MSSMA, MSSMB and MSSMC are given, then the specified parameters are used to 
calculate all the masses and decay modes with the ISASUSY package assuming the minimal 
supersymmetric extension of the standard model (MSSM). There are no default values, so 
you must specify values for each MSSMi, i=A-C. MSSMD can optionally be used to set the 
second generation squark and slepton parameters; if it is omitted, then the first generation 
ones are used. MSSME can optionally be used to set the U(l) and SU(2) gaugino masses; 
if it is omitted, then the grand unification values are used. The parameters and the use of 
the MSSM is preserved if the title is SAME. FORCE can be used to override the calculated 
branching ratios. 

The MSSM option also generates charginos and neutralinos with cross sections based on 
the MSSM mixing angles in addition to squarks and sleptons. These can be selected with 
JETTYPEi; the complete list of supersymmetric options is: 

'GLSS' , 

'UPSSL' , 'UBSSL' , 'DNSSL' , 'DBSSL' , 'STSSL' , 'SBSSL' , 'CHSSL' , 'CBSSL' , 
'BTSS1' , 'BBSS1' , 'TPSS1' , 'TBSS1' , 



44 



'UPSSR' , 'UBSSR' , 'DNSSR' , 'DBSSR' , 'STSSR' , 'SBSSR' , 'CHSSR' , 'CBSSR' , 
'BTSS2' , 'BBSS2 ' , 'TPSS2 ' , 'TBSS2' , 

'W1SS+' , ' W1SS- ' , 'W2SS+' , 'W2SS-' , 'Z1SS' , 'Z2SS' , 'Z3SS' , 'Z4SS' , 

'NUEL' , ' ANUEL ' , >EL-> , >EL+' , ' NUML ' , ' ANUML ' ,MUL-' , 'MUL+' , ' NUTL ' , 

' ANUTL ' , 'TAU1- ' , 'TAU1+' , ' ER- ' , 'ER+' , 'MUR-' , ' MUR+ ' , 'TAU2-' , 'TAU2+' , 

'ZO' , 'HLO' , 'HHO' , 'HAO' , 'H+' , 'H-' , 

' SQUARKS ' , ' GAUGINOS ' , ' SLEPTONS ' , ' ALL ' . 

Note that mixing between L and R stop states results in 1 (light) and 2 (heavy) stop, sbottom 
and stau eigenstates, which depend on the input parameters of left- and right- scalar masses, 
plus A terms, fi and tan (3. The last four JETTYPE's generate respectively all allowed 
combinations of squarks and antisquarks, all combinations of charginos and neutralinos, all 
combinations of sleptons and sneutrinos, and all SUSY particles. 

For SUSY Higgs pair production or associated production in E+E-, select the appropriate 
JETTYPE's, e.g. 

JETTYPE1 
'ZOV 
JETTYPE2 
'HLO'/ 

As usual, this gives only half the cross section. For single production of neutral SUSY Higgs 
in pp and pp reactions, use the HIGGS process together with the MSSMi, SUGRA, or GMSB 
keywords. You must specify one and only one Higgs type using 

HTYPE 

'HLO' or 'HHO' or 'HAO'/ ««< One only! 

If no QMH range is given, one is calculated using M ± 5T for the selected Higgs. Decays 
into quarks, leptons, gauge bosons, lighter Higgs bosons, and SUSY particles are generated 
using the on-shell branching ratios from ISASUSY. You can use JETTYPEi to select the 
allowed Higgs modes and WMODEi to select the allowed decays of W and Z bosons. Since 
heavy SUSY Higgs bosons couple weakly to W pairs, WW fusion and WW scattering are 
not included. 

SUGRA can be used instead of MSSMi to generate MSSM decays with parameters de- 
termined from m , miw, A , tan/3, and sgn/i = ±1 in the minimal supergravity framework. 
The NUSUGi keywords can optionally be used to specify additional parameters for non- 
universal SUGRA models. Similarly, the GMSB keyword is used to specify the A, M m , A 5 , 
tan (3, sgn/i = ±1, and C grav parameters of the minimal Gauge Mediated SUSY Breaking 
model. GMSB2 can optionally be used to specify additional parameters of non-minimal 
GMSB models. 

WHIGGS is used to generate W plus neutral Higgs events. For the Standard Model 
the JETTYPE is HIGGS. If any of the SUSY models is specified, then the appropriate SUSY 
Higgs type should be used, most likely HLO. In either case WMODEi is used to specify the W 
decay modes. The Higgs is treated as a particle; its decay modes can be set using FORCE. 

The FORCE keyword requires special care. Its list must contain the numerical particle 
IDENT codes, e.g. 
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FORCE 

140,130,-120/ 

The charge-conjugate mode is also forced for its antiparticle. Thus the above example forces 
both D° — > K + n~ and D° — > fT~7r + . If only a specific decay is wanted one should use the 
FORCE1 command; e.g. 

F0RCE1 

140,130,-120/ 

only forces D° — > K + n~ . 

To force a heavy quark decay one must generally separately force each hadron containing 
it. If the decay is into three leptons or quarks, then the real or virtual W propagator is 
inserted automatically. Since Version 7.30, top and fourth generation quarks are treated as 
particles and decayed directly rather than first being made into hadrons. Thus for example 

F0RCE1 
6,-12,11,5/ 

forces all top quarks to decay into an positron, neutrino and a b-quark (which will be 
hadronized). For the physical top mass, the positron and neutrino will come from a real W. 
Note that forcing t — > W + b and W + — > e + v e does not give the same result; the first uses the 
correct V — A matrix element, while the second decays the W according to phase space. 

Forced modes included in the decay table will automatically be put into the correct order. 
Modes not listed in the decay table are allowed, but caution is advised because a wrong decay 
mode can cause an infinite loop or other unexpected effects. 

FORCE (FORCE1) can be called at most 20 (40) times in any run plus all subsequent 
'SAME' runs. If it is called more than once for a given parent, all calls are listed, and the 
last call is used. Note that FORCE applies to particles only, but that for gamma, W+, W-, 
Z0 and supersymmetric particles the same IDENT codes are used both as jet types and as 
particles. 

The default parton distributions are fit CTEQ3L from the CTEQ Collaboration using 
lowest order QCD. The CTEQ and the older EHLQ and Duke-Owens distributions can be 
selected using the STRUC keyword. 

If PDFLIB support is enabled (see Section 4), then any of the distributions in the PDFLIB 
compilation by H. Plothow-Besch can be selected using the PDFLIB keyword and giving the 
proper parameters, which are identical to those described in the PDFLIB manual and are 
simply passed to the routine PDFSET. For example, to select fit 29 (CTEQ3L) by the CTEQ 
group, leaving all other parameters with their default values, use 

PDFLIB 

' CTEQ ' , 29D0/ 

Note that the fit-number and the other parameters are of type DOUBLE PRECISION 
(REAL on 64-bit machines). There is no internal passing of parameters except for those 
which control the printing of messages. 
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7 Output 



The output tape or file contains three types of records. A beginning record is written by 
a call to ISAWBG before generating a set of events; an event record is written by a call to 
ISAWEV for each event; and an end record is written for each run by a call to ISAWND. 
These subroutines load the common blocks described below into a single 

COMMON/ZEVEL/ZEVEL ( 1024) 

and write it out when it is full. A subroutine RDTAPE, described in the next section, inverts 
this process so that the user can analyze the event. 

ZEVEL is written out to TAPEj by a call to BUFOUT. For the CDC version IF = 
PAIRPAK is selected; BUFOUT first packs two words from ZEVEL into one word in 

C0MM0N/ZV0UT/ZV0UT(512) 

using subroutine PAIRPAK and then does a buffer out of ZVOUT to TAPEj. Typically at 
least two records are written per event. For all other computers IF=STDIO is selected, and 
ZEVEL is written out with a standard FORTRAN unformatted write. 



7.1 Beginning Record 

At the start of each run ISAWBG is called. It writes out the following common blocks: 

COMMON/DYLIM/QMIN , QMAX , QTMIN, QTMAX , YWMIN , YWMAX , XWMIN , XWMAX , THWMIN , 

2 THWMAX,PHWMIN,PHWMAX 

3 ,SETLMQ(12) 
SAVE /DYLIM/ 
LOGICAL SETLMQ 
EQUIVALENCE (BLIMl(l) ,QMIN) 

REAL QMIN , QMAX , QTMIN , QTMAX , YWMIN , YWMAX , XWMIN , XWMAX , THWMIN , 

+ THWMAX,PHWMIN,PHWMAX,BLIM1(12) 

QMIN, QMAX = W mass limits 

QTMIN,QTMAX = W q t limits 

YWMIN, YWMAX = W rj rapidity limits 

XWMIN,XWMAX = W x F limits 

THWMIN, THWMAX = W 6 limits 

PHWMIN,PHWMAX = W <p limits 

COMMON/ IDRUN/ IDVER, IDG (2) , IEVT, IEVGEN 
SAVE /IDRUN/ 

INTEGER IDVER , IDG , IEVT , IEVGEN 
IDVER = program version 

IDG(l) = run date (10000 xmonth+ 100 xday+year) 
IDG(2) = run time (10000 xhour+lOOxminute+second) 
IEVT = event number 



47 



COMMON/ JETLIM/PMIN (3) ,PMAX(3) ,PTMIN(3) ,PTMAX(3) ,YJMIN(3) ,YJMAX(3) 

1 , PHIMIN (3) , PHIMAX (3) , X JMIN (3) , X JMAX (3) , THMIN (3) , THMAX (3) 

2 ,SETLMJ(36) 
SAVE /JETLIM/ 

EQUIVALENCE (BLIMS(l) ,PMIN(1)) 
LOGICAL SETLMJ 

C0MM0N/FIXPAR/FIXP(3) ,FIXPT(3) ,FIXYJ(3) ,FIXPHI(3) ,FIXXJ(3) 

2 , FIXQM , FIXQT , FIXYW , FIXXW , FIXPHW 
LOGICAL FIXQM, FIXQT , FIXYW , FIXXW , FIXPHW 
LOGICAL FIXP , FIXPT , FIXY J , FIXPHI , FIXXJ 
C0MM0N/SGNPAR/CTHS(2,3) ,THS(2,3) ,YJS(2,3) ,XJS(2,3) 
REAL PMIN , PMAX , PTMIN , PTMAX , Y JMIN , Y JMAX , PHIMIN , PHIMAX , X JMIN , 

+ X JMAX , THMIN , THMAX , BLIMS (36) , CTHS , THS , Y JS , X JS 

PMIN,PMAX = jet momentum limits 

PTMIN, PTMAX = jet p t limits 

Y JMIN ,Y JMAX = jet r) rapidity limits 

PHIMIN,PHIMAX = jet <p limits 

THMIN, THMAX = jet 9 limits 

COMMON/KEYS/ IKEYS , KEYON , KEYS (10) 

COMMON/XKEYS/REAC 

SAVE /KEYS/ , /XKEYS/ 

LOGICAL KEYS 

LOGICAL KEYON 

CHARACTER*8 REAC 

INTEGER IKEYS 

KEYON = normally TRUE, FALSE if no good reaction 
KEYS = TRUE if reaction I is chosen 

1 for TWOJET 

2 for E+E- 

3 for DRELLYAN 

4 for MINBIAS 

5 for SUPERSYM 

6 for WPAIR 

REAC = character reaction code 



COMMON/PRIMAR/N JET , SCM , HALFE , ECM , IDIN (2) , NEVENT , NTRIES , NSIGMA 
SAVE /PRIMAR/ 

INTEGER N JET , IDIN , NEVENT , NTRIES , NSIGMA 
REAL SCM, HALFE, ECM 
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NJET = number of jets per event 

SCM = square of com energy 

HALFE = beam energy 

ECM = com energy 

IDIN = ident code for initial beams 

NEVENT = number of events to be generated 

NTRIES = maximum number of tries for good jet parameters 

NSIGMA = number of extra events to determine SIGF 

INTEGER MXGOQ 
PARAMETER (MXG0Q=85) 

C0MM0N/qiQ2/G0q(MXG0q,3) ,G0ALL(3) ,G0DY(4) ,STDDY,G0WW(25,2) , 
$ALLWW(2) ,G0WM0D(25,3) 
SAVE /Q1Q2/ 

LOGICAL GOQ , GOALL , GODY , STDDY , GOWW , ALLWW , GOWMOD 



GOQ(I,K) 



GOALL(K) 
GODY(I) 
1= 1 2 3 4 
GM W+ W- 
STDDY 
GOWW(I,K) 
ALLWW(K) 



TRUE if quark type I allowed for jet k 
1=1234 5_6 789 10 11 12 13 

=^-guuddssccbbti 
I = 14 15 16 17 18 19 20 21 22 23 24 25 

=>• v e v e e" e + fT fi + v T v T r~ r" 

TRUE if all jet types allowed 
TRUE if W type I is allowed. 



Z0 



REAL 

ALAM 
ALAM2 
CUTJET 
ISTRUC 



= TRUE if standard DRELLYAN 
= TRUE if I is allowed in the decay of K for WPAIR. 
= TRUE if all allowed in the decay of K for WPAIR. 

COMMON/qCDPAR/ ALAM , ALAM2 , CUTJET , ISTRUC 
SAVE /qCDPAR/ 
INTEGER ISTRUC 

ALAM, ALAM2, CUT JET 

QCD scale A 
QCD scale A 2 

cutoff for generating secondary partons 
3 for Eichten (EHLQ), 
= 4 for Duke (DO) 
= 5 for CTEQ 2L 
= 6 for CTEQ 3L 
= -999 for PDFLIB 

COMMON/qLMASS/ AMLEP ( 100) , NqLEP , NMES , NBARY 
SAVE /qLMASS/ 

I NTEGER NqLEP , NMES , NBARY 
REAL AMLEP 



AMLEP(6:£ 



t,y,x masses, only elements written 
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7.2 Event Record 

For each event ISAWEV is called. It writes out the following common blocks: 

COMMON/FINAL/ NKINF , SIGF , ALUM , ACCEPT , NRECS 

SAVE /FINAL/ 

INTEGER NKINF, NRECS 

REAL SIGF, ALUM, ACCEPT 

SIGF = integrated cross section, only element written 

COMMON/ IDRUN/IDVER, IDG (2) , IEVT, IEVGEN 
SAVE /IDRUN/ 

INTEGER IDVER , IDG , IEVT , IEVGEN 

IDVER = program version 
IDG = run identification 
IEVT = event number 

COMMON/ JETPAR/P (3) ,PT(3) ,YJ(3) ,PHI(3) ,XJ(3) ,TH(3) , CTH (3) , STH(3) 

1 , JETTYP (3) ,SHAT,THAT,UHAT,QSQ,X1,X2,PBEAM(2) 

2 , QMW , QW , QTW , YW , XW , THW , QTMW , PHI W , SHAT 1 , THAT 1 , UHAT 1 , JWTYP 

3 ,ALFQSQ,CTHW,STHW,QOW 

4 ,INITYP(2) ,ISIGS,PBEAMS(5) 
SAVE /JETPAR/ 

INTEGER JETTYP , JWTYP , INITYP , ISIGS 



REAL P, PT,YJ, PHI, XJ,TH, CTH, STH, SHAT, THAT, UHAT, QSQ, XI, X2, 

+ PBEAM , QMW , QW , QTW , YW , XW , THW , QTMW , PHIW , SHAT1 , THAT1 , UHAT1 , 

+ ALFQSQ , CTHW , STHW , QOW , PBEAMS 

P = jet momentum \p\ 

PT = jet pt 

YJ = jet 77 rapidity 

PHI = jet 

XJ = jet xf 

TH = jet 9 

CTH = jet cos(#) 

STH = jet sin(0) 

JETTYP = jet type. The code is listed under /Q1Q2/ above 
continued. . . 
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PITT Am TITT Am T TTT Am 

SHA1,1HA1,UHA1 


— 


hard scattering s, t, u 


QSQ 


— 


effective Q z 






initial parton rr^ 


PBEAM 


— — 


1 j 

remaining beam momentum 






TTT" 

W mass 


ow 




\Af mompntnm 

V V LLL\JLLL\jLL\j L1X11 


OTW 




T/l/ T"Tfl TT QA/'PTQP TT10TTlPTTt"lTm 
K v ul dllo VC1 DC llJAJlllCllu U.111 


YW 

-L V V 




V V L CXiYJ IUIIj y 


xw 

-/V V V 




W Tzr 


THW 




W 9 


(al 1 IVl VV 




1 '2 1 /02 
V ?t,W + V 


PHIW 


= 


TV 


OUT A T1 TTT A rni t ttt a rp i 

SHA1 1, THAI 1,UHA1 1 




invariants for TT decay 


JWTYP 




W type. The code is listed under /Q1Q2/ above. 


ALFQSQ 




QCD coupling a s (Q 2 ) 


CTHW 




W cos(#) 


STHW 




W sin(0) 


QOW 




W energy 



INTEGER MXJSET,JPACK 

PARAMETER (MXJSET=400 , JPACK=1000) 

C0MM0N/JETSET/NJSET,PJSET(5,MXJSET) , JORIG(MXJSET) , JTYPE(MXJSET) , 
$JDCAY(MXJSET) 
SAVE /JETSET/ 

INTEGER NJSET , JORIG , JTYPE , JDCAY 
REAL PJSET 

NJSET = number of partons 

PJSET(1,I) = p x of parton I 
PJSET(2,I) = p y of parton I 
PJSET(3,I) = p z of parton I 
PJSET(4,I) = p of parton I 
PJSET(5,I) = mass of parton I 

JORIG(I) = JPACK*JET+K if I is a decay product of K. 

IF K=0 then I is a primary parton. 

(JET = 1,2,3 for final jets.) 

(JET = 11,12 for initial jets.) 
JTYPE(I) = IDENT code for parton I 

JDCAY(I) = JPACK*K1+K2 if Kl and K2 are decay products of I. 

If JDCAY(I)=0 then I is a final parton 
MXJSET = dimension for /JETSET/ arrays. 
JPACK = packing integer for /JETSET/ arrays. 

INTEGER MXSIGS , IOPAK 

PARAMETER (MXSIGS=3000 , I0PAK=100) 
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COMMON/JETSIG/SIGMA,SIGS(MXSIGS) ,NSIGS , INOUT(MXSIGS) , SIGEVT 

SAVE /JETSIG/ 

INTEGER NSIGS,INOUT 

REAL SIGMA, SIGS, SIGEVT 

SIGMA = cross section summed over types 

SIGS(I) = cross section for reaction I (not written) 

NSIGS = number of nonzero cross sections (not written) 

INOUT(I) = packed partons for process I (not written) 

MXSIGS = dimension for JETSIG arrays (not written) 

SIGEVT = partial cross section for selected channel 

INTEGER MXPTCL , IPACK 

PARAMETER (MXPTCL=4000 , IPACK=10000) 

COMMON/PARTCL/NPTCL , PPTCL ( 5 , MXPTCL ) , IORIG (MXPTCL) , IDENT (MXPTCL) 
l.IDCAY (MXPTCL) 
SAVE /PARTCL/ 

INTEGER NPTCL, IORIG, IDENT, IDCAY 
REAL PPTCL 



NPTCL 

PPTCL(1,I 

PPTCL(2,I 

PPTCL(3,I 

PPTCL(4,I 

PPTCL(5,I 

IORIG(I) 



IDENT(I) 
IDCAY(I) 

MXPTCL 
IPACK 



number of particles 
p x for particle I 
p y for particle I 
p z for particle I 
Po for particle I 
mass for particle I 

IPACK* JET+K if I is a decay product of K. 

- (IPACK* JET+K) if I is a primary particle from 

parton K in /JETSET/. 

if I is a primary beam particle. 

(JET = 1,2,3 for final jets.) 

(JET = 11,12 for initial jets.) 

IDENT code for particle I 

IPACK*K1+K2 if decay products are K1-K2 inclusive. 
If IDCAY(I)=0 then particle I is stable, 
dimension for /PARTCL/ arrays, 
packing integer for /PARTCL/ arrays. 



C0MM0N/PINITS/PINITS(5,2) ,IDINIT(2) 
SAVE /PINITS/ 
INTEGER IDINIT 
REAL PINITS 
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PINITS(1,I 
PINITS(2,I 
PINITS(3,I 
PINITS(4,I 
PINITS(5,I 
IDINIT(I) 



p x for initial parton I 
p y for initial parton I 
p z for initial parton I 
p for initial parton I 
mass for initial parton I 
IDENT for initial parton I 



INTEGER MXJETS 
PARAMETER (MXJETS=10) 

COMMON/PJETS/PJETS (5, MXJETS) , ID JETS (MXJETS) , QW JET (5) , IDENTW 
$,PPAIR(5,4) ,IDPAIR(4) , JPAIR(4) , NPAIR, IFRAME (MXJETS) 
SAVE /PJETS/ 



INTEGER 


ID JETS , IDENTW , IDPAIR , JPAIR , NPAIR , IFRAME 


REAL 


PJETS , QW JET , PPAIR 


PJETS(1,I) = 


for iet I 


PJETS(2,I) = 


j0„ for iet I 

fy iWJ - 


PJETS(3,I) = 


v? for iet I 


PJETS(4,I) = 


Po for jet I 


PJETS(5,I) = 


mass for jet I 


IDJETS(I) = 


IDENT code for jet I 


QWJET(l) = 


p x for W 


QWJET(2) = 


p y for W 


QWJET(3) = 


p z for W 


QWJET(4) = 


Pq for W 


QWJET(5) = 


mass for W 


IDENTW 


IDENT CODE for W 


PPAIR(1,I) = 


p x for WPAIR decay product I 


PPAIR(2,I) = 


p y for WPAIR decay product I 


PPAIR(3,I) = 


p z for WPAIR decay product I 


PPAIR(4,I) = 


Po for WPAIR decay product I 


PPAIR(5,I) = 


mass for WPAIR decay product I 


IDPAIR(I) = 


IDENT code for WPAIR product I 


JPAIR(I) 


JETTYPE code for WPAIR product I 


NPAIR 


2 for VT /± 7 events, 4 for WW events 



COMMON/TOTALS/NKINPT , NWGEN , NKEEP , SUMWT , WT 
SAVE /TOTALS/ 

INTEGER NKINPT, NWGEN, NKEEP 
REAL SUMWT, WT 

NKINPT = number of kinematic points generated. 

NWGEN = number of W+jet events accepted. 

NKEEP = number of events kept. 

SUMWT = sum of weighted cross sections. 

WT = current weight. (SIGMAxWT = event weight.) 
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COMMON/WSIG/SIGLLQ 
SAVE /WSIG/ 
REAL SIGLLQ 

SIGLLQ = cross section for W decay. 

Of course irrelevant common blocks such as /WSIG/ for TWO JET events are not written 
out. 



7.3 End Record 

At the end of a set ISAWND is called. It writes out the following common block: 

COMMON/FINAL/NKINF , SIGF , ALUM , ACCEPT , NRECS 

SAVE /FINAL/ 

INTEGER NKINF, NRECS 

REAL S I GF , ALUM , ACCEPT 

NKINF = number of points generated to calculate SIGF 
SIGF = integrated cross section for this run 

ALUM = equivalent luminosity for this run 
ACCEPT = ratio of events kept over events generated 
NRECS = number of physical records for this run 

Events within a given run have uniform weight. Separate runs can be combined together 
using SIGF/NEVENT as the weight per event. This gives a true cross section in mb units. 

The user can replace subroutines ISAWBG, ISAWEV, and ISAWND to write out the 
events in a different format or to update histograms using HBOOK or any similar package. 
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8 File Reading 



The FORTRAN instruction 

CALL RDTAPE (IDEV,IFL) 

will read a beginning record, an end record or an event (which can be more than one record). 
IDEV is the tape number and 

IFL=0 for a good read, 
IFL=-1 for an end of file. 

The information is restored to the common blocks described above. The type of record is 
contained in 

COMMON/RECTP/IRECTP , IREC 

SAVE /RECTP/ 

INTEGER IRECTP , IREC 

IRECTP = 100 for an event record 
IRECTP = 200 for a beginning record 
IRECTP = 300 for an end record 

IREC = no. of physical records in event record, otherwise 

The parton momenta from the primary hard scattering are contained in /PJETS/. The 
parton momenta generated by the QCD cascade are contained in /JETSET/. The hadron 
momenta both from the QCD jets and from the beam jets are contained in /PARTCL/. 
The final hadron momenta and the associated pointers should be used to calculate the jet 
momenta, since they are changed both by the QCD cascade and by hadronization. Particles 
with IDCAY=0 are stable, while the others are resonances. 

The weight per event needed to produce a weighted histogram in millibarn units is 
SIGF/NEVENT. The integrated cross section SIGF is calculated by Monte Carlo integra- 
tion during the run for the given kinematic limits and JETTYPE, WTYPE, and WMODE 
selections. Any of three methods can be used to find the value of SIGF: 

(1) The current value, which is written out with each event, can be used. To prevent 
enormous fluctuations at the beginning of a run, NSIGMA extra primary parton events are 
generated first. The default value, NSIGMA = 20, gives negligible overhead but may not be 
large enough for good accuracy. 

(2) The value SIGF calculated with the full statistics of the run can be obtained by 
reading through the tape until an end record (IRECTP=300) is found. After SIGF is saved 
with a different name, the first event record for the run can be found by backspacing the 
tape NRECS times. 

(3) Unweighted histograms can be made for the run and the weight added after the end 
record is found. An implementation of this using special features of HBO OK is contained in 
ISAPLT. 

The functions AMASS(IDENT), CHARGE(IDENT), and LABEL(IDENT) are available 
to determine the mass, charge, and character label in A8 format. Subroutine FLAVOR 
returns the quark content of any hadron and may be useful to convert IDENT codes to 
other schemes. CALL PRTEVT(O) prints an event. 
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9 Decay Table 



ISAJET uses an external table of decay modes. Particles can be put into the table in 
arbitrary order, but all modes for each particle must be grouped together. The table is 
rewound and read in before each run with a READ* format. Each entry must have the form 

IDENT , ITYPE , CBR , ID1 , ID2 , ID3 , ID4 , ID5/ 

where IDENT is the code for the parent particle, ITYPE specifies the type of decay - 
ITYPE = 1 for strong, ITYPE = 2 for electromagnetic, ITYPE = 3 for weak — CBR is the 
cumulative branching ratio, and ID1,. . . ,ID5 are the IDENT codes for the decay products. 
The parent IDENT must be positive; the charge conjugate mode is used for the antiparticle. 
The values of CBR must of course be positive and monotonically increasing for each mode, 
with the last value being 1.00 for each parent IDENT. The last parent IDENT code must be 
zero. Care should be taken in adding new modes, since there is no checking for validity. In 
some cases order is important; note in particular that quarks and gluons must always appear 
last so that they can be removed and fragmented into hadrons. 

After the decay table there is a space for text, which is printed under the heading of 
ISAJET NOTES for each run. It must be terminated by a line containing ////. This 
feature is used at BNL to announce new versions and other changes. 

The decay table is contained in patch ISADECAY. 
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10 IDENT Codes 



ISAJET uses a numerical ident code for particle types. Quarks and leptons are numbered 
in order of mass: 

UP =1 NUE = 11 

DN = 2 E- = 12 

ST =3 NUM = 13 

CH = 4 MU- = 14 

BT =5 NUT = 15 

TP =6 TAU- = 16 

with a negative sign for antiparticles. Arbitrary conventions are: 

GL =9 
GM = 10 
KS = 20 
KL =-20 
W+ = 80 
Z0 = 90 

The supersymmetric particle IDENT codes distinguish between the partners of left and right 
handed fermions and include the Higgs sector of the minimal supersymmetric model: 

UPSSL . . . TPSS1 =21 ... 26 

NUEL . . . TAU1- =31 ... 36 

UPSSR . . . TPSS2 =41 ... 46 

NUER . . . TAU2- =51 ... 56 
GLSS = 29 

Z1SS = 30 Z2SS = 40 

Z3SS = 50 Z4SS = 60 

W1SS+ = 39 W2SS+ = 49 



HL0 = 82 HH0 = 83 

HA0 =84 H+ = 86 

The code for a meson is a compound integer +-JKL, where J.LE.K are the quarks and L is 
the spin. The sign is for the J quark. Glueball IDENT codes have not been selected, but 
the choice GL=9 clearly allows 990, 9990, etc. Flavor singlet mesons are ordered by mass, 

PI0 = 110 
ETA = 220 
ETAP = 330 
ETAC = 440 

which is natural for the heavy quarks. Similarly, the code for a baryon is a compound integer 
+-IJKL formed from the three quarks I,J,K and a spin label L=0,1. The code for a diquark 
is +-IJ00. Additional states are distinguished by a fifth integer, e.g., 
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A1+ = 10121 



These and a few J=2 mesons are used in some of the B decays. 

A routine PRTLST is provided to print out a complete list of valid IDENT codes and 
associated information. The usage is CALL PRTLST(LUN, AMY, AMX) where LUN is the 
unit number and AMY and AMX are the masses of the Y and X quarks respectively. This 
routine should be linked with the ISAJET library and with ALDATA. 

The complete list of ident codes follows. (Hadrons containing t quarks are defined but 
are no longer listed since the t quark is treated as a particle.) 



IDENT 


LABEL 




MASS 


CHARGE 


1 


UP 




.30000E+00 


.67 


-1 


UB 




.30000E+00 


-.67 


2 


DN 




.30000E+00 


-.33 


-2 


DB 




.30000E+00 


.33 


3 


ST 




.50000E+00 


-.33 


-3 


SB 




.50000E+00 


.33 


4 


CH 




. 16000E+01 


.67 


-4 


CB 




. 16000E+01 


-.67 


5 


BT 




.49000E+01 


-.33 


-5 


BB 




.49000E+01 


.33 


6 


TP 




. 17500E+03 


.67 


-6 


TB 




. 17500E+03 


-.67 


9 


GL 


0, 




0.00 


10 


GM 


0. 




0.00 


11 


NUE 


0. 




0.00 


-11 


ANUE 


0, 




0.00 


12 


E- 




.51100E-03 


-1.00 


-12 


E+ 




.51100E-03 


1.00 


13 


NUM 


0. 




0.00 


-13 


ANUM 


0, 




0.00 


14 


MU- 




. 10566E+00 


-1.00 


-14 


MU+ 




. 10566E+00 


1.00 


15 


NUT 


0. 




0.00 


-15 


ANUT 


0. 




0.00 


16 


TAU- 




. 18070E+01 


-1.00 


-16 


TAU+ 




. 18070E+01 


1.00 


20 


KS 




.49767E+00 


0.00 


-20 


KL 




.49767E+00 


0.00 


21 


UPSSL 


none 


0.67 
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-21 UBSSL none -0.67 

22 DNSSL none -0.33 
-22 DBSSL none 0.33 

23 STSSL none -0.33 

23 SBSSL none 0.33 

24 CHSSL none 0.67 
-24 CBSSL none -0.67 

25 BTSS1 none -0.33 
-25 BBSS1 none 0.33 

26 TPSS1 none 0.67 
-26 TBSS1 none -0.67 

29 GLSS none 0.00 

30 Z1SS none 0.00 

31 NUEL none 0.00 
-31 ANUEL none 0.00 

32 EL- none -1.00 
-32 EL+ none +1.00 

33 NUML none 0.00 
-33 ANUML none 0.00 

34 MUL- none -1.00 
-34 MUL+ none +1.00 

35 NUTL none 0.00 
-35 ANUTL none 0.00 

36 TAU1- none -1.00 
-36 TAU1+ none -1.00 

39 W1SS+ none 1.00 
-39 W1SS- none -1.00 

40 Z2SS none 0.00 

41 UPSSR none 0.67 
-41 UBSSR none -0.67 

42 DNSSR none -0.33 
-42 DBSSR none 0.33 

43 STSSR none -0.33 

43 SBSSR none 0.33 

44 CHSSR none 0.67 
-44 CBSSR none -0.67 

45 BTSS2 none -0.33 
-45 BBSS2 none 0.33 

46 TPSS2 none 0.67 
-46 TBSS2 none -0.67 
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49 


W2SS+ 


none 


1. 


.00 


-49 


W2SS- 


none 


-1. 


.00 


50 


Z3SS 


none 


0, 


.00 


51 


NUER 


none 


0, 


.00 


-51 


ANUER 


none 


0, 


.00 


52 


ER- 


none 


-1, 


.00 


-52 


ER+ 


none 


+1, 


.00 


53 


NUMR 


none 


0, 


.00 


-53 


ANUMR 


none 


0, 


.00 


54 


MUR- 


none 


-1 , 


.00 


-54 


MUR+ 


none 


+1, 


.00 


55 


NUTR 


none 


0. 


.00 


-55 


ANUTR 


none 


0. 


.00 


56 


TAU2- 


none 


-1. 


.00 


-56 


TAU2+ 


none 


-1. 


.00 
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Z4SS 


none 


0, 


.00 


80 


W+ 


.80200E+02 


1, 


.00 


81 


HIGGS 


.80200E+02 


0. 


.00 


82 


HLO 


none 


0, 


.00 


83 


HHO 


none 


0, 


.00 


84 


HAO 


none 


0, 


.00 


86 


H+ 


none 


1, 


.00 


90 


ZO 


.91190E+02 


0, 


.00 


91 


GVSS 





0, 


.00 


110 


PIO 


. 13496E+00 


0, 


.00 


120 


PI+ 


. 13957E+00 


1, 


.00 


120 


PI- 


. 13957E+00 


-1, 


.00 


220 


ETA 


. 54745E+00 


0, 


.00 


130 


K+ 


.49367E+00 


1, 


.00 


130 


K- 


.49367E+00 


-1. 


.00 


230 


KO 


.49767E+00 


0. 


.00 


230 


AKO 


.49767E+00 


0. 


.00 


330 


ETAP 


.95760E+00 


0. 


.00 


140 


ADO 


.18645E+01 


0. 


.00 


140 


DO 


.18645E+01 


0. 


.00 


240 


D- 


.18693E+01 


-1, 


.00 


240 


D+ 


.18693E+01 


1, 


.00 


340 


F- 


.19688E+01 


-1, 


.00 


340 


F+ 


.19688E+01 


1, 


.00 



440 


ETAC 


.29788E+01 


0, 


.00 


150 


UB. 


.51700E+01 


1. 


.00 


-150 


BU. 


.51700E+01 


-1. 


.00 


250 


DB. 


.51700E+01 


0, 


.00 


-250 


BD. 


.51700E+01 


. 


. 00 


350 


SB. 


.53700E+01 


0, 


.00 


-350 


BS. 


.53700E+01 


0, 


.00 


450 


CB. 


.64700E+01 


1, 


.00 


-450 


BC. 


.64700E+01 


-1, 


.00 


550 


BB. 


.97700E+01 


0, 


.00 


111 


RH00 


.76810E+00 


0, 


.00 


121 


RHO+ 


.76810E+00 


1, 


.00 


-121 


RHO- 


.76810E+00 


-1. 


.00 


221 


OMEG 


.78195E+00 


0. 


.00 


131 


K*+ 


.89159E+00 


1. 


.00 


-131 


K*- 


.89159E+00 


-1. 


.00 


231 


K*0 


.89610E+00 


0. 


.00 


-231 


AK*0 


.89610E+00 





.00 


331 


PHI 


.10194E+01 


0. 


.00 


141 


AD*0 


.20071E+01 


0, 


.00 


-141 


D*0 


.20071E+01 


0, 


.00 


241 


D*- 


.20101E+01 


-1, 


.00 


-241 


D*+ 


.20101E+01 


1, 


.00 


341 


F*- 


.21103E+01 


-1, 


.00 


-341 


F*+ 


.21103E+01 


1, 


.00 


441 


JPSI 


.30969E+01 


0, 


.00 


151 


UB* 


.52100E+01 


1, 


.00 


-151 


BU* 


.52100E+01 


-1, 


.00 


251 


DB* 


.52100E+01 


0, 


.00 


-251 


BD* 


.52100E+01 


. 


. 00 


351 


SB* 


.54100E+01 


0, 


.00 


-351 


BS* 


.54100E+01 


0. 


.00 


451 


CB* 


.65100E+01 


1, 


.00 


-451 


BC* 


.65100E+01 


-1. 


.00 


551 


UPSL 


.98100E+01 


0. 


.00 


112 


F2 


.12750E+01 


0. 


.00 


132 


K2*+ 


.14254E+01 


1. 


.00 


-132 


K2*- 


.14254E+01 


-1, 


.00 


232 


K2*0 


.14324E+01 


0, 


.00 


-232 


AK2*0 


.14324E+01 


0, 


.00 


.0110 


FO 


.98000E+00 


0. 


.00 
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10111 


A10 


.12300E+01 


0, 


.00 


10121 


A1+ 


.12300E+01 


1, 


.00 


-10121 


Al- 


.12300E+01 


-1, 


.00 


10131 


K1+ 


.12730E+01 


1, 


.00 


-10131 


Kl- 


. 12730E+01 


-1 . 


.00 


10231 


K10 


.12730E+01 


0, 


.00 


-10231 


AK10 


. 12730E+01 


0, 


.00 


30131 


Kl*+ 


.14120E+01 


1, 


.00 


-30131 


Kl*- 


.14120E+01 


-1. 


.00 


30231 


K1*0 


.14120E+01 


0, 


.00 


-30231 


AK1*0 


.14120E+01 


0, 


.00 


10441 


PSK2S) 


.36860E+01 


0. 


.00 


20440 


CHIO 


.34151E+01 


0. 


.00 


20441 


CHI1 


.35105E+01 


0. 


.00 


20442 


CHI2 


.35662E+01 


0. 


.00 


1120 


P 


.93828E+00 


1, 


.00 


-1120 


AP 


.93828E+00 


-1, 


.00 


1220 


N 


.93957E+00 


0. 


.00 


-1220 


AN 


.93957E+00 


0, 


.00 


1130 


S+ 


.11894E+01 


1, 


.00 


-1130 


AS- 


.11894E+01 


-1, 


.00 


1230 


SO 


.11925E+01 


0, 


.00 


-1230 


ASO 


.11925E+01 


0, 


.00 


2130 


L 


.11156E+01 


0, 


.00 


-2130 


AL 


.11156E+01 


0, 


.00 


2230 


s- 


.11974E+01 


-1, 


.00 


-2230 


AS+ 


.11974E+01 


1, 


.00 


1330 


XIO 


.13149E+01 


0, 


.00 


-1330 


AXIO 


.13149E+01 


0, 


.00 


2330 


XI- 


.13213E+01 


-1. 


.00 


-2330 


AXI+ 


.13213E+01 


1. 


.00 


1140 


SC++ 


.24527E+01 


2. 


.00 


-1140 


ASC— 


.24527E+01 


-2. 


.00 


1240 


SC+ 


.24529E+01 


1. 


.00 


-1240 


ASC- 


.24529E+01 


-1, 


.00 


2140 


LC+ 


.22849E+01 


1, 


.00 


-2140 


ALC- 


.22849E+01 


-1, 


.00 


2240 


SCO 


.24525E+01 


0, 


.00 


-2240 


ASCO 


.24525E+01 


0. 


.00 
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1340 


use. 


.25000E+01 


1. 


.00 


-1340 


AUSC. 


.25000E+01 


-1. 


.00 


3140 


sue. 


.24000E+01 


1, 


.00 


-3140 


ASUC. 


.24000E+01 


-1, 


.00 


2340 


DSC. 


.25000E+01 


0, 


.00 


-2340 


ADSC. 


.25000E+01 


0, 


.00 


3240 


SDC. 


.24000E+01 


0, 


.00 


-3240 


ASDC. 


.24000E+01 


0, 


.00 


3340 


SSC. 


.26000E+01 


0, 


.00 


-3340 


ASSC. 


.26000E+01 


0. 


.00 


1440 


UCC. 


.35500E+01 


2, 


.00 


-1440 


AUCC. 


.35500E+01 


-2. 


.00 


2440 


DCC. 


.35500E+01 


1. 


.00 


-2440 


ADCC. 


.35500E+01 


-1. 


.00 


3440 


sec. 


.37000E+01 


1. 


.00 


-3440 


ASCC. 


.37000E+01 


-1. 


.00 


1150 


UUB. 


.54700E+01 


1. 


.00 


-1150 


AUUB. 


.54700E+01 


-1. 


.00 


1250 


UDB. 


.54700E+01 


0, 


.00 


-1250 


AUDB. 


.54700E+01 


0, 


.00 


2150 


DUB. 


.54700E+01 


0, 


.00 


-2150 


ADUB. 


.54700E+01 


0, 


.00 


2250 


DDB. 


.54700E+01 


-1, 


.00 


-2250 


ADDB. 


.54700E+01 


1, 


.00 


1350 


USB. 


.56700E+01 


0, 


.00 


-1350 


AUSB. 


.56700E+01 


0, 


.00 


3150 


SUB. 


.56700E+01 


0, 


.00 


-3150 


ASUB. 


.56700E+01 


0, 


.00 


2350 


DSB. 


.56700E+01 


-1. 


.00 


-2350 


ADSB. 


.56700E+01 


1, 


.00 


3250 


SDB. 


.56700E+01 


-1, 


.00 


-3250 


ASDB. 


.56700E+01 


1, 


.00 


3350 


SSB. 


.58700E+01 


-1, 


.00 


-3350 


ASSB. 


.58700E+01 


1, 


.00 


1450 


UCB. 


.67700E+01 


1. 


.00 


-1450 


AUCB. 


.67700E+01 


-1, 


.00 


4150 


CUB. 


.67700E+01 


1. 


.00 


-4150 


ACUB. 


.67700E+01 


-1. 


.00 


2450 


DCB. 


.67700E+01 


0. 


.00 


-2450 


ADCB. 


.67700E+01 


0. 


.00 


4250 


CDB. 


.67700E+01 


0, 


.00 


-4250 


ACDB. 


.67700E+01 


0, 


.00 


3450 


SCB. 


.69700E+01 


0, 


.00 


-3450 


ASCB. 


.69700E+01 


0, 


.00 
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4350 


CSB. 


.69700E+01 


0, 


.00 


-4350 


ACSB. 


.69700E+01 


0, 


.00 


4450 


CCB. 


.80700E+01 


1. 


.00 


-4450 


ACCB. 


.80700E+01 


-1, 


.00 


1550 


UBB. 


. 10070E+02 


0, 


.00 


-1550 


AUBB. 


. 10070E+02 


0, 


.00 


2550 


DBB. 


. 10070E+02 


-1 , 


.00 


-2550 


ADBB . 


. 10070E+02 


1, 


.00 


3550 


SBB . 


. 10270E+02 


-1, 


.00 


-3550 


ASBB. 


. 10270E+02 


1, 


.00 


4550 


CBB. 


.11370E+02 


0, 


.00 


-4550 


ACBB . 


.11370E+02 


0, 


.00 


1111 


DL++ 


.12320E+01 


2. 


.00 


-1111 


ADL— 


.12320E+01 


-2. 


.00 


1121 


DL+ 


.12320E+01 


1. 


.00 


-1121 


ADL- 


.12320E+01 


-1. 


.00 


1221 


DLO 


.12320E+01 


0. 


.00 


-1221 


ADLO 


.12320E+01 





.00 


2221 


DL- 


.12320E+01 


-1, 


.00 


-2221 


ADL+ 


.12320E+01 


1, 


.00 


1131 


S*+ 


.13823E+01 


1, 


.00 


-1131 


AS*- 


.13823E+01 


-1, 


.00 


1231 


S*0 


.13820E+01 


0, 


.00 


-1231 


AS*0 


.13820E+01 


0. 


.00 


2231 


s*- 


.13875E+01 


-1, 


.00 


-2231 


AS*+ 


.13875E+01 


1, 


.00 


1331 


XI*0 


.15318E+01 


0, 


.00 


-1331 


AXI*0 


.15318E+01 


0, 


.00 


2331 


XI*- 


.15350E+01 


-1. 


.00 


-2331 


AXI*+ 


.15350E+01 


1, 


.00 


3331 


OM- 


.16722E+01 


-1, 


.00 


-3331 


AOM+ 


.16722E+01 


1, 


.00 


1141 


UUC* 


.26300E+01 


2, 


.00 


-1141 


AUUC* 


.26300E+01 


-2. 


.00 


1241 


UDC* 


.26300E+01 


1. 


.00 


-1241 


AUDC* 


.26300E+01 


-1. 


.00 


2241 


DDC* 


.26300E+01 


0. 


.00 


-2241 


ADDC* 


.26300E+01 


0. 


.00 


1341 


use* 


.27000E+01 


1. 


.00 


-1341 


AUSC* 


.27000E+01 


-1, 


.00 


2341 


DSC* 


.27000E+01 


0, 


.00 


-2341 


ADSC* 


.27000E+01 


0, 


.00 


3341 


SSC* 


.28000E+01 


0, 


.00 
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-3341 


ASSC* 


.28000E+01 


0, 


.00 


1441 


UCC* 


.37500E+01 


2, 


.00 


-1441 


AUCC* 


.37500E+01 


-2, 


.00 


2441 


DCC* 


.37500E+01 


1. 


.00 


-2441 


ADCC* 


.37500E+01 


-1. 


.00 


3441 


sec* 


.39000E+01 


1, 


.00 


-3441 


ASCC* 


.39000E+01 


-1, 


.00 


4441 


CCC* 


.48000E+01 


2, 


.00 


-4441 


ACCC* 


.48000E+01 


-2, 


.00 


1151 


UUB* 


.55100E+01 


1, 


.00 


-1151 


AUUB* 


.55100E+01 


-1, 


.00 


1251 


UDB* 


.55100E+01 


0, 


.00 


-1251 


AUDB* 


.55100E+01 


0, 


.00 


2251 


DDB* 


.55100E+01 


-1. 


.00 


-2251 


ADDB* 


.55100E+01 


1. 


.00 


1351 


USB* 


.57100E+01 


0. 


.00 


-1351 


AUSB* 


.57100E+01 


0. 


.00 


2351 


DSB* 


.57100E+01 


-1. 


.00 


-2351 


ADSB* 


.57100E+01 


1. 


.00 


3351 


SSB* 


.59100E+01 


-1, 


.00 


-3351 


ASSB* 


.59100E+01 


1, 


.00 


1451 


UCB* 


.68100E+01 


1, 


.00 


-1451 


AUCB* 


.68100E+01 


-1, 


.00 


2451 


DCB* 


.68100E+01 


0, 


.00 


-2451 


ADCB* 


.68100E+01 


0, 


.00 


3451 


SCB* 


.70100E+01 


0, 


.00 


-3451 


ASCB* 


.70100E+01 


0. 


.00 


4451 


CCB* 


.81100E+01 


1, 


.00 


-4451 


ACCB* 


.81100E+01 


-1, 


.00 


1551 


UBB* 


.10110E+02 


0, 


.00 


-1551 


AUBB* 


.10110E+02 


0, 


.00 


2551 


DBB* 


.10110E+02 


-1, 


.00 


-2551 


ADBB* 


.10110E+02 


1, 


.00 


3551 


SBB* 


.10310E+02 


-1, 


.00 


-3551 


ASBB* 


.10310E+02 


1. 


.00 


4551 


CBB* 


.11410E+02 


0. 


.00 


-4551 


ACBB* 


.11410E+02 


0. 


.00 


5551 


BBB* 


.14710E+02 


-1. 


.00 


-5551 


ABBB* 


.14710E+02 


1, 


.00 


1100 


UUO. 


.60000E+00 





.67 


-1100 


AUUO. 


.60000E+00 


-o. 


.67 


1200 


UDO. 


.60000E+00 


0. 


.33 
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-1200 


AUDO. 


.60000E+00 


-0, 


.33 


2200 


DDO . 


.60000E+00 


-0. 


.67 


-2200 


ADDO. 


.60000E+00 


0, 


.67 


1300 


USO. 


.80000E+00 


0, 


.33 


-1300 


AUSO. 


.80000E+00 


-0, 


.33 


2300 


DSO. 


.80000E+00 


-0. 


.67 


-2300 


ADSO. 


.80000E+00 


0, 


.67 


3300 


SSO. 


.10000E+01 


-0. 


.67 


-3300 


ASSO. 


.10000E+01 


0, 


.67 


1400 


UCO. 


.19000E+01 


1. 


.33 


-1400 


AUCO. 


.19000E+01 


-1. 


.33 


2400 


DCO. 


.19000E+01 


0, 


.33 


-2400 


ADCO. 


.19000E+01 


-o. 


.33 


3400 


SCO. 


.21000E+01 


0. 


.33 


-3400 


ASCO. 


.21000E+01 


-o. 


.33 


4400 


CCO. 


.32000E+01 


1. 


.33 


-4400 


ACCO. 


.32000E+01 


-1. 


.33 


1500 


UBO. 


.49000E+01 


0. 


.33 


-1500 


AUBO. 


.49000E+01 


-o. 


.33 


2500 


DBO. 


.49000E+01 


-o. 


.67 


-2500 


ADBO. 


.49000E+01 


0. 


.67 


3500 


SBO. 


.51000E+01 


-o. 


.67 


-3500 


ASBO. 


.51000E+01 


0. 


.67 


4500 


CBO. 


.65000E+01 


0, 


.33 


-4500 


ACBO. 


.65000E+01 


-o. 


.33 


5500 


BBO. 


.98000E+01 


-o. 


.67 


-5500 


ABBO. 


.98000E+01 


0, 


.67 
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11 Higher Order Processes 



Higher order processes can be generated either by the QCD evolution or by supplying partons 
from an external generator. 

Frequently it is interesting to generate higher-order processes with a particular branching 
in the QCD evolution or with a particular particle or group of particles being produced from 
the fragmentation. Examples include 

1. Branching of jets into heavy quarks (e.g., g — > b + b); 

2. Decay of such a heavy quark into a lepton or neutrino; 

3. Radiation of a photon, W, or Z from a jet. 

It is important to realize that all of the cross sections and the QCD evolution in ISAJET 
are based on leading-log QCD, so generating such processes does not give the correct higher 
order QCD cross sections or "K factors" , even though it may produce better agreement with 
them in some cases. 

ISAJET does produce events with particular topologies which in many cases are the most 
important effect of higher order processes. In the heavy quark example, the lowest order 
process 

g + g^Q + Q 

produces back-to-back heavy quark pairs, whereas the splitting process 

g + g^g + g, g^Q + Q 

produces collinear pairs. Such collinear pairs are essential to obtain agreement with experi- 
mental data on bb production, and they often are the dominant background for processes of 
interest. 

Branchings such as the emission of a heavy quark pair, a photon, or a W ± or Z° are 
rare, and since they may occur at any step in the evolution, one cannot force them to occur. 
Therefore, generation of such events is very slow. M. Delia Negra (UAl) suggested first 
doing rt\ QCD evolutions for each hard scattering and rejecting events without the desired 
partons, then doing n 2 fragmentations for each successful evolution. This generates the 
equivalent of nin 2 events for each hard scattering, so the cross section must be divided by 
7ii7i 2 . This algorithm can speed up the generation of g — > b + b splitting by a factor of ten 
for n\ = n 2 = 10. 

Since the evolution and fragmentation steps are executed nin 2 times even if good events 
are found, a single hard scattering can lead to multiple events. This does not change the 
inclusive cross sections, but it does mean that the fluctuations may be larger than expected. 
Hence it is important to choose the numbers ri\ and n 2 carefully. 

The following entities are used in ISAJET for generating events with multiple evolution 
and fragmentation: 

NEVENT: The number of primary hard scatterings to be generated. Set as usual on the 
input line with the energy. 
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SIGF: The cross section for the selected hard scatterings divided by ri\ x n 2 . Hence the 
correct weight is SIGF/NEVENT, just as for normal running. (The cross section printed at 
the end of a run does not contain this factor.) 

NEVOLVE: The number n\ of evolutions per hard scattering. This should never be set 
unless you supply a REJJET function. Do not confuse this with NOEVOLVE. 

NHADRON: The number n 2 of fragmentations for a given evolution. This should never be 
set unless you supply a REJFRG function. Do not confuse this with NOHADRON. 

REJJET: A logical function which if true causes the evolution to be rejected. The user 
must supply one to make the selections which he wants. The default always .FALSE, but 
includes an example comment. 

REJFRG: A logical function which if true causes the fragmentation to be rejected. The 
user must supply one to make the selections which he wants. The default always .FALSE, 
but includes an example as a comment. 

Note that one can also use function EDIT to make a final selection of the events. Of course 
ISAJET must be relinked if EDIT, REJJET or REJFRG is modified. 

At the end of a run, the jet cross section, the cross section for the selected events, and the 
number and fraction of events selected are printed. The cross section SIGF stored internally 
is divided by rt\ x n 2 so that if the events are used to make histograms, then the correct 
weight per event is 

SIGF/NEVENT 

just as for normal events. Of course NEVENT now has a different meaning; it is in gen- 
eral larger than the number of events in the file but might be smaller if NEVOLVE and 
NHADRON are badly chosen. 

NEVOLVE and NHADRON are set as parameters in the input. One wants to choose 
them to give better acceptance of the primary hard scatterings but not to give multiple 
events for one hard scattering. For lepton production from heavy quarks the values 

NEVOLVE 
10/ 

NHADRON 
10/ 

seem appropriate, giving reasonable efficiency For radiation of photons from jets, NEVOLVE 
can be somewhat larger but NHADRON should be one, and REJFRG should always return 
.FALSE., since the selection is just on the parton process, not on the hadronization. 

The loops over evolutions and fragmentations are done inside of subroutine ISAEVT 
and are always executed the same number of times even though ISAEVT returns after each 
generated event. Logical flag OK signals a good event, and logical flag DONE signals that 
the run is finished. If you control the event generation loop yourself, you should make use 
of these flags as in the following extract from subroutine ISAJET: 

IL00P=0 
101 CONTINUE 

IL00P=IL00P+1 
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CALL ISAEVT(ILOOP,OK,DONE) 
IF (OK) CALL ISAWEV 
IF (.NOT. DONE) GO TO 101 

Otherwise you may get the wrong weights. 

It is possible to supply to ISAJET events with partons generated by some other program 
that may have more accurate matrix elements for higher order processes. Because any such 
calculation must involve cutoffs ISAJET assumes that the partons were generated imposing 
some R cutoff, where R = \/(j) 2 + r/ 2 , and some E t cutoff. Given that information ISAJET 
will generate initial state radiation partons only below the Et cutoff and final state radiation 
inside the R cutoff. The external partons can be supplied to ISAJET by calls to 2 subroutines. 
To initialize ISAJET for externally supplied partons, use 

CALL INISAP (CMSE , REACTION , BEAMS , WZ , NDCAYS , DCAYS , ETMIN , RCONE , OK) 

where the inputs are 
CMSE = center of mass energy 

REACTION = reaction (only TWOJET and DRELLYAN are 
implemented so far) 

BEAMS(2) = chose 'P ' or 'AP' 

ETMIN = minimum ET of supplied partons 

RCONE = minimum cone (R) between supplied partons 

WZ = option 'W, 'Z', or ' ' no Ws or Z's 

NDCAYS = number of decay options (if 0, assume decay has 
already been done) 

DCAYS = list of particles W or Z can decay into 

and the output is 

OK = TRUE if initialization is possible 
Then for each event use 

CALL IPARTNS (NPRTNS , IDS , PRTNS , IDQ , WEIGHT , WZDK) 

where the inputs are 
NPRTNS = number of partons, < 10 

IDS(NPRTNS) = ids of final partons 

PRTNS (4,NPRTNS) = parton 4 vectors 
IDQ(2) = ids of initial partons 

WEIGHT = weight 

WZDK = if true last 2 partons are from W,Z decay 

Further QCD radiation is then generated consistent with ETMIN and RCONE, and the 
partons are fragmented into hadrons as usual. If RCONE is set to a value greater than 1.5 
no cone restriction is applied during parton evolution. 



69 



12 ISASUSY: Decay Modes in the Minimal Supersym- 
metric Model 



The code in patch ISASUSY of ISAJET calculates decay modes of supersymmetric parti- 
cles based on the work of H. Baer, M. Bisset, M. Drees, D. Dzialo (Karatas), X. Tata, J. 
Woodside, and their collaborators. The calculations assume the minimal supersymmetric 
extension of the standard model. The user specifies the gluino mass, the pseudoscalar Higgs 
mass, the Higgsino mass parameter /i, tan/5, the soft breaking masses for the first and third 
generation left-handed squark and slepton doublets and right-handed singlets, and the third 
generation mixing parameters A t , A b , and A T . Supersymmetric grand unification is assumed 
by default in the chargino and neutralino mass matrices, although the user can optionally 
specify arbitrary U(l) and £77(2) gaugino masses at the weak scale. The first and second 
generations are assumed by default to be degenerate, but the user can optionally specify 
different values. These inputs are then used to calculate the mass eigenstates, mixings, and 
decay modes. 

Most calculations are done at the tree level, but one-loop results for gluino loop decays, 
H — > 77 and H — > gg, loop corrections to the Higgs mass spectrum and couplings, and 
leading-log QCD corrections to H — > qq are included. The Higgs masses have been calculated 
using the effective potential approximation including both top and bottom Yukawa and 
mixing effects. Mike Bisset and Xerxes Tata have contributed the Higgs mass, couplings, and 
decay routines. Manuel Drees has calculated several of the three-body decays including the 
full Yukawa contribution, which is important for large tan(beta). Note that e+e- annihilation 
to SUSY particles and SUSY Higgs bosons have been included in ISAJET versions > 7.11. 
ISAJET versions > 7.22 include the large tan (3 solution as well as non-degenerate sfermion 
masses. 

Other processes may be added in future versions as the physics interest warrants. Note 
that the details of the masses and the decay modes can be quite sensitive to choices of 
standard model parameters such as the QCD coupling ALFA3 and the quark masses. To 
change these, you must modify subroutine SSMSSM. By default, ALFA3=.12. 

All the mass spectrum and branching ratio calculations in ISASUSY are performed by a 
call to subroutine SSMSSM. Effective with version 7.23, the calling sequence is 

SUBROUTINE SSMSSM (XMG , XMU , XMHA , XTANB , XMQ 1 , XMDR , XMUR , 
$XML 1 , XMER , XMQ2 , XMSR , XMCR , XML2 , XMMR ,XMQ3, XMBR , XMTR , 
$XML3 , XMLR, XAT , XAB , XAL , XM1 , XM2 , XMT , I ALLOW) 

where the following are taken to be independent parameters: 
XMG = gluino mass 
XMU = n = SUSY Higgs mass 



XMHA 
XTANB 



—2 * mi of Baer et al. 
pseudo-scalar Higgs mass 
tan/3, ratio of vev's 
1/R (of old Baer- Tata notation). 
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XMQ1 = qi soft mass, 1st generation 

XMDR = d r mass, 1st generation 

XMUR = u r mass, 1st generation 

XML1 = i\ mass, 1st generation 

XMER = e r mass, 1st generation 

XMQ2 = qi soft mass, 2nd generation 

XMSR = s r mass, 2nd generation 

XMCR = c r mass, 2nd generation 

XML2 = £i mass, 2nd generation 

XMMR = jl r mass, 2nd generation 

XMQ3 = qi soft mass, 3rd generation 

XMBR = b r mass, 3rd generation 

XMTR = t r mass, 3rd generation 

XML3 = l\ mass, 3rd generation 

XMTR = f r mass, 3rd generation 

XAT = stop trilinear term A t 

XAB = sbottom trilinear term A^ 

XAL = stau trilinear term A T 

XM1 = U(l) gaugino mass 

= computed from XMG if I 1E19 

XM2 = SU(2) gaugino mass 

= computed from XMG if I 1E19 

XMT = top quark mass 
The variable IALLOW is returned: 

IALLOW = 1 if Z1SS is not LSP, otherwise 
All variables are of type REAL except IALLOW, which is INTEGER, and all masses are in 
GeV. The notation is taken to correspond to that of Haber and Kane, although the Tata 
Lagrangian is used internally. All other standard model parameters are hard wired in this 
subroutine; they are not obtained from the rest of ISAJET. The theoretically favored range 
of these parameters is 

50 < M(g) < 2000 GeV 
50 < M(q) < 2000 GeV 
50 < M{1) < 2000 GeV 
-1000 < /i < 1000 GeV 
1 < tan/5 < m t /m b 
M(t) « 175 GeV 
50 < M(A) < 2000 GeV 
M(t,),M(t r )<M(g) 
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M(b r ) ~ M(q) 
-1000 < A t < 1000 GeV 
-1000 < A b < 1000 GeV 

It is assumed that the lightest supersymmetric particle is the lightest neutralino Zi, the 
lighter stau fi, or the gravitino G in GMSB models. Some choices of the above parameters 
may violate this assumption, yielding a light chargino or light stop squark lighter than Z\. 
In such cases SSMSSM does not compute any branching ratios and returns IALLOW = I. 

SSMSSM does not check the parameters or resulting masses against existing experimental 
data. SSTEST provides a minimal test. This routine is called after SSMSSM by ISAJET 
and ISASUSY and prints suitable warning messages. 

SSMSSM first calculates the other SUSY masses and mixings and puts them in the 
common block /SSPAR/: 



c 


SUSY parameters 






c 


AMGLSS 


_ 


gluino mass 


c 


AMULSS 

1 11 1 \J J— 1 K-J 





wo— left sauark mass 


c 


AMELSS 





lef t-selectron mass 


c 


AMERSS 

ill ij—ix kj 


= 


ricrht - sleDton mass 


c 


AMNiSS 





sneutrino mass for generation i 


c 


TW0M1 




Hi ctctq i n n masR = — mn 

11 J. c. CLO 111U U1CLO O 111 LI 


c 


RV2V1 





ratio v2/vl of vev's 


c 


AMTLSS , AMTRSS 





left, right stop masses 


c 


AMT1SS,AMT2SS 




light, heavy stop masses 


c 


AMBLSS , AMBRSS 




left, right sbottom masses 


c 


AMB1SS,AMB2SS 




light, heavy sbottom masses 


c 


AMLLSS , AMLRSS 




left, right stau masses 


c 


AML1SS,AML2SS 




light, heavy stau masses 


c 


AMZiSS 




signed mass of Zi 


c 


ZMIXSS 




Zi mixing matrix 


c 


AMWiSS 




signed Wi mass 


c 


GAMMAL , GAMMAR 




Wi left, right mixing angles 


c 


AMHL , AMHH , AMHA 




neutral Higgs hO, HO, AO masses 


c 


AMHC 




charged Higgs H+ mass 


c 


ALFAH 




Higgs mixing angle 


c 


AAT 




stop trilinear term 


c 


THETAT 




stop mixing angle 


c 


AAB 




sbottom trilinear term 


c 


THETAB 




sbottom mixing angle 


c 


AAL 




stau trilinear term 


c 


THETAL 




stau mixing angle 


c 


AMGVSS 




gravitino mass 



C0MM0N/SSPAR/ AMGLSS , AMULSS , AMURSS , AMDLSS , AMDRSS , AMSLSS 
$ , AMSRSS , AMCLSS , AMCRSS , AMBLSS , AMBRSS , AMB1SS , AMB2SS 
$ , AMTLSS , AMTRSS , AMT1SS , AMT2SS , AMELSS , AMERSS , AMMLSS , AMMRSS 
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$ , AMLLSS , AMLRSS , AML1SS , AML2SS , AMN1SS , AMN2SS , AMN3SS 
$ , TW0M1 , RV2V1 , AMZ1SS , AMZ2SS , AMZ3SS , AMZ4SS , ZMIXSS (4 ,4) 
$,AMW1SS,AMW2SS 

$ , GAMMAL , GAMMAR , AMHL , AMHH , AMHA , AMHC , ALFAH , AAT , THETAT 
$ , AAB , THETAB , AAL , THETAL , AMGVSS 

REAL AMGLSS , AMULSS , AMURSS , AMDLSS , AMDRSS , AMSLSS 
$ , AMSRSS , AMCLSS , AMCRSS , AMBLSS , AMBRSS , AMB1SS , AMB2SS 
$ , AMTLSS , AMTRSS , AMT1SS , AMT2SS , AMELSS , AMERSS , AMMLSS , AMMRSS 
$ , AMLLSS , AMLRSS , AML1SS , AML2SS , AMN1SS , AMN2SS , AMN3SS 
$ , TWOM1 , RV2V1 , AMZ1SS , AMZ2SS , AMZ3SS , AMZ4SS , ZMIXSS 
$,AMW1SS,AMW2SS 

$ , GAMMAL , GAMMAR , AMHL , AMHH , AMHA , AMHC , ALFAH , AAT , THETAT 
$ , AAB , THETAB , AAL , THETAL , AMGVSS 
REAL AMZISS(4) 

EQUIVALENCE (AMZISS(l) ,AMZ1SS) 
SAVE /SSPAR/ 

It then calculates the widths and branching ratios and puts them in the common block 
/SSMODE/: 



C MXSS = maximum number of modes 

C NSSMOD = number of modes 

C ISSMOD = initial particle 

C JSSMOD = final particles 

C GSSMOD = width 

C BSSMOD = branching ratio 



INTEGER MXSS 
PARAMETER (MXSS=1000) 

COMMON/SSMODE/NSSMOD , ISSMOD (MXSS) , JSSMOD (5 , MXSS) , GSSMOD (MXSS) 
$, BSSMOD (MXSS) 
INTEGER NSSMOD, ISSMOD, JSSMOD 
REAL GSSMOD, BSSMOD 
SAVE /SSMODE/ 

Decay modes for a given particle are not necessarily adjacent in this common block. Note 
that the branching ratio calculations use the full matrix elements, which in general will give 
nonuniform distributions in phase space, but this information is not saved in /SSMODE/. 
In particular, the decays H-^Z + Z*^Z + f + f give no indication that the // mass is 
strongly peaked near the upper limit. 

All IDENT codes are defined by parameter statements in the PATCHY keep sequence 
SSTYPE: 

C SM ident code definitions. These are standard ISAJET but 

C can be changed. 

INTEGER IDUP , IDDN , IDST , IDCH , IDBT , IDTP 

INTEGER IDNE , IDE , IDNM , IDMU , IDNT , IDTAU 



73 



INTEGER IDGL , IDGM , IDW , IDZ 

PARAMETER (IDUP=1 , IDDN=2 , IDST=3 , IDCH=4 , IDBT=5 , IDTP=6) 

PARAMETER (IDNE=1 1 , IDE=12 , IDNM=13 , IDMU=14 , IDNT=15 , IDTAU=16) 

PARAMETER (IDGL=9 , IDGM=10 , IDW=80 , IDZ=90) 
C SUSY ident code definitions. They are chosen to be similar 

C to those in versions < 6.50 but may be changed. 

INTEGER ISUPL , ISDNL , ISSTL , ISCHL , ISBT1 , ISTP1 

I NTEGER I SNEL , I SEL , I SNML , I SMUL , ISNTL, ISTAU1 

INTEGER ISUPR, ISDNR, ISSTR, ISCHR, ISBT2, ISTP2 

INTEGER ISNER, ISER, ISNMR, ISMUR, ISNTR, ISTAU2 

INTEGER ISZ1 , ISZ2 , ISZ3 , ISZ4 , ISW1 , ISW2 , ISGL 

INTEGER ISHL,ISHH,ISHA,ISHC 

INTEGER ISGRAV 

PARAMETER (ISUPL=21 , ISDNL=22 , ISSTL=23 , ISCHL=24 , ISBT1=25 , ISTP1=26) 
PARAMETER (ISNEL=31 , ISEL=32 , ISNML=33 , ISMUL=34 , ISNTL=35 , ISTAU1=36) 
PARAMETER (ISUPR=41 , ISDNR=42 , ISSTR=43 , ISCHR=44 , ISBT2=45 , ISTP2=46) 
PARAMETER (ISNER=51 , ISER=52 , ISNMR=53 , ISMUR=54 , ISNTR=55 , ISTAU2=56) 
PARAMETER (ISGL=29) 

PARAMETER (ISZ1=30 , ISZ2=40 , ISZ3=50 , ISZ4=60 , ISW1=39 , ISW2=49) 
PARAMETER (ISHL=82 , ISHH=83 , ISHA=84, ISHC=86) 
PARAMETER (ISGRAV=91) 

These are based on standard ISAJET but can be changed to interface with other generators. 
Since masses except the t mass are hard wired, one should check the kinematics for any 
decay before using it with possibly different masses. 

Instead of specifying all the SUSY parameters at the electroweak scale using the MSSMi 
commands, one can instead use the SUGRA parameter to specify in the minimal supergravity 
framework the common scalar mass m , the common gaugino mass m 1( / 2 , and the soft trilinear 
SUSY breaking parameter A at the GUT scale, the ratio tan/3 of Higgs vacuum expectation 
values at the electroweak scale, and sgn //, the sign of the Higgsino mass term. The NUSUGi 
keywords allow one to break the assumption of universality in various ways. NUSUGI sets the 
gaugino masses; NUSUG2 sets the A terms; NUSUG3 sets the Higgs masses; NUSUG4 sets the 
first generation squark and slepton masses; and NUSUG5 sets the third generation masses. 

The renormalization group equations are solved iteratively using Runge-Kutta numerical 
integration to determine the weak scale parameters from the GUT scale ones: 

1. The RGE's are run from the weak scale Mz up to the GUT scale, where a± = «2, 
taking all thresholds into account. We use two loop RGE equations for the gauge 
couplings only. 

2. The GUT scale boundary conditions are imposed, and the RGE's are run back to Mz, 
again taking thresholds into account. 

3. The masses of the SUSY particles and the values of the soft breaking parameters B 
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and mu needed for radiative symmetry are computed, e.g. 



, , . Ml - Ml tan 2 8 9 , 

"< Mz >= V-i M|/2 



These couplings are frozen out at the scale y M(t L )M(t R ). 

4. The 1-loop radiative corrections are computed. 

5. The process is then iterated until stable results are obtained. 

This is essentially identical to the procedure used by several other groups. Other possible 
constraints such as b-tau unification and limits on proton decay have not been included. 

An alternative to the SUGRA model is the Gauge Mediated SUSY Breaking (GMSB) 
model of Dine and Nelson, Phys. Rev. D48, 1277 (1973); Dine, Nelson, Nir, and Shirman, 
Phys. Rev. D53, 2658 (1996). In this model SUSY is broken dynamically and communicated 
to the MSSM through messenger fields at a messenger mass scale M m much less than the 
Planck scale. If the messenger fields are in complete representations of SU(5), then the 
unification of couplings suggested by the LEP data is preserved. The simplest model has 
a single 5 + 5 messenger sector with a mass M m and and a SUSY-breaking VEV F rn of its 
auxiliary field F . Gauginos get masses from one-loop graphs proportional to A m = F m /M m 
times the appropriate gauge coupling af, sfermions get squared- masses from two- loop graphs 
proportional to A m times the square of the appropriate ctj. If there are N 5 messenger fields, 
the gaugino masses and sfermion masses- squared each contain a factor of N 5 . 

The parameters of the GMSB model implemented in ISAJET are 

• A m = F m /M m : the scale of SUSY breaking, typically 10-100 TeV; 

• M m > A m : the messenger mass scale, at which the boundary conditions for the renor- 
malization group equations are imposed; 

• N 5 : the equivalent number of 5 + 5 messenger fields. 

• tan/3: the ratio of Higgs vacuum expectation values at the electroweak scale; 
sgn/i = ±1: the sign of the Higgsino mass term; 



• Cgrav > 1: the ratio of the gravitino mass to the value it would have had if the only 
SUSY breaking scale were F m . 

The solution of the renormalization group equations is essentially the same as for SUGRA; 
only the boundary conditions are changed. In particular it is assumed that electroweak 
symmetry is broken radiatively by the top Yukawa coupling. 

In GMSB models the lightest SUSY particle is always the nearly massless gravitino G. 
The phenomenology depends on the nature of the next lightest SUSY particle (NLSP) and 
on its lifetime to decay to a gravitino. The NLSP can be either a neutralino Xi or a slepton 
T\ . Its lifetime depends on the gravitino mass, which is determined by the scale of SUSY 
breaking not just in the messenger sector but also in any other hidden sector. If this is set 
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by the messenger scale F m , i.e., if C gr3iV ~ 1, then this lifetime is generally short. However, 
if the messenger SUSY breaking scale F m is related by a small coupling constant to a much 
larger SUSY breaking scale F b , then C grav 3> 1 and the NLSP can be long-lived. The correct 
scale is not known, so C grav should be treated as an arbitrary parameter. More complicated 
GMSB models may be run by using the GMSB2 keyword. 

Patch ISASSRUN of ISAJET provides a main program SSRUN and some utility programs 
to produce human readable output. These utilities must be rewritten if the IDENT codes in 
/SSTYPE/ are modified. To create the stand-alone version of ISASUSY with SSRUN, run 
YPATCHY on isajet.car with the following cradle (with & replaced by +): 

&USE,* ISASUSY. Select all code 

&USE,NOCERN. No CERN Library 

&USE,IMPNONE. Use IMPLICIT NONE 

&EXE. Write everything to ASM 

&PAM,T=C. Read PAM file 

&QUIT. Quit 

Compile, link, and run the resulting program, and follow the prompts for input. Patch 
ISASSRUN also contains a main program SUGRUN that reads the minimal SUGRA, non- 
universal SUGRA, or GMSB parameters, solves the renormalization group equations, and 
calculates the masses and branching ratios. To create the stand-alone version of ISASUGRA, 
run YPATCHY with the following cradle: 

&USE, * ISASUGRA . Select all code 

&USE,NOCERN. No CERN Library 

&USE,IMPNONE. Use IMPLICIT NONE 

&EXE. Write everything to ASM 

&PAM. Read PAM file 

&QUIT. Quit 

The documentation for ISASUSY and ISASUGRA is included with that for ISAJET. 

ISASUSY is written in ANSI standard Fortran 77 except that IMPLICIT NONE is used 
if +USE,IMPNONE is selected in the Patchy cradle. All variables are explicitly typed, and 
variables starting with I,J,K,L,M,N are not necessarily integers. All external names such as 
the names of subroutines and common blocks start with the letters SS. Most calculations 
are done in double precision. If +USE,NOCERN is selected in the Patchy cradle, then the 
Cernlib routines EISRS1 and its auxiliaries to calculate the eigenvalues of a real symmetric 
matrix and DDILOG to calculate the dilogarithm function are included. Hence it is not 
necessary to link with Cernlib. 

The physics assumptions and details of incorporating the Minimal Supersymmetric Model 
into ISAJET have appeared in a conference proceedings entitled "Simulating Supersymmetry 
with ISAJET 7.0/ISASUSY 1.0" by H. Baer, F. Paige, S. Protopopescu and X. Tata; this 
has appeared in the proceedings of the workshop on Physics at Current Accelerators and 
Supercolliders, ed. J. Hewett, A. White and D. Zeppenfeld, (Argonne National Laboratory, 
1993). Detailed references may be found therein. Users wishing to cite an appropriate source 
may cite the above report. 
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13 Changes in Recent Versions 



This section contains a record of changes in recently released versions of ISAJET, taken 
from the memoranda distributed to users. Note that the released version numbers are not 
necessarily consecutive. 

13.1 Version 7.40, October 1998 

A new process WHIGGS generates W ± + H and Z + H events for both the Standard Model 
and SUSY models and also Higgs pair production for SUSY models. The types and W 
decay modes are selected with JETTYPE and WMODE as for WPAIR events. This process 
is of particular interest for producing fairly light Higgs bosons at the Tevatron. See the 
documentation for more details. 

Some non-minimal GMSB models can be generated using a new keyword GMSB2. The 
optional parameters are an extra factor between the gaugino and scalar masses, shifts in the 
Higgs masses, a D-term proportional to hypercharge, and independent numbers of messenger 
fields for the three gauge groups. The documentation gives more details and references. 

The default for SUGRA models has been changed to use a s (Mz) = 0.118, the experimen- 
tal value. This means that the couplings do not exactly unify at the GUT scale, presumably 
because of the effects of heavy particles. The keyword AL3UNI can be used to select exact 
unification, which produces too large a value for a s (M z ). 

A number of three-body slepton decays that occur through left-right mixing are now 
included. These are obviously small but might compete with gravitino decays. In particular, 
a decay like JIr — > f\vv might lead to a wrong momentum measurement in the muon system. 
So far we have found no case in which this is probable. 

The new release also includes a separate Unix tar file mcpp.tar containing C++ code to 
read a standard ISAJET output file and copy all the information into C++ classes. The tar 
file contains makefiles for Software Release Tools, documentation, and examples as well as 
the code. 

13.2 Version 7.37, April 1998 

Version 7.37 incorporates Gauge Mediated SUSY Breaking models for the first time. In 
these models, SUSY is broken in a hidden sector at a relatively low scale, and the masses 
of the MSSM fields are then produced through ordinary gauge interactions with messenger 
fields. The parameters of the GMSB model in ISAJET are M m , the messenger mass scale; 
A m = F m /M m , where F m is the SUSY breaking scale in the messenger sector; N 5 , the 
number of messenger fields; the usual tan/3 and sgn/x; and C grav > 1, a factor which scales 
the gravitino mass and hence the lifetime for the lightest MSSM particle to decay into it. 

GMSB models have a light gravitino G as the lightest SUSY particle. The phenomenology 
of the model depends mainly on the nature of the next lightest SUSY particle, a x? or a fi, 
which changes with the number N 5 of messengers. The phenomenology also depends on the 
lifetime for the Xi ~^ 01 ^1 ~~ * Gt decay; this lifetime can be short or very long. All the 
relevant decays are included except for jl — > vufi, which is very suppressed. 
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The keyword MGVTNO allows the user to independently input a gravitino gravitino 
mass for the MSSM option. This allows studies of SUGRA (or other types) of models where 
the gravitino is the LSP. 

Version 7.37 also contains an extension of the SUGRA model with a variety of non- 
universal gaugino and sfermion masses and A terms at the GUT scale. This makes it possible 
to study for example, how well the SUGRA assumptions can be tested. 

Two significant bugs have also been corrected. The decay modes for B* mesons were 
missing from the decay table since Version 7.29 and have been restored. A sign error in the 
interference term for chargino production has been corrected, leading to a larger chargino 
pair cross section at the Tevatron. 

13.3 Version 7.32, November 1997 

This version makes several corrections in various chargino and neutralino widths, thus chang- 
ing the branching ratios for large tan/3. For x%, f° r example, the Xi^b branching ratio is 
decreased significantly, while the Xi T+T ~ one is increased. Thus the SUGRA phenomenology 
for tan (3 ~ 30 is modified substantially. 

The new version also fixes a few bugs, including a possible numerical precision problem 
in the Drell-Yan process at high mass and qr- It also includes a missing routine for the Zebra 
interface. 

13.4 Version 7.31, August 1997 

Version fixes a couple of bugs in Version 7.29. In particular, the JETTYPE selection did not 
work correctly for supersymmetric Higgs bosons, and there was an error in the interactive 
interface for MSSM input. Since these could lead to incorrect results, users should replace 
the old version. We thank Art Kreymer for finding these problems. 

Since top quarks decay before they have time to hadronize, they are now put directly 
onto the particle list. Top hadrons (tu, td, etc.) no longer appear, and FORCE should be 
used directly for the top quark, i.e. 

FORCE 

6,11,-12,5/ 

The documentation has been converted to LaTeX. Run either LaTeX 2.09 or LaTeX 2e 
three times to resolve all the forward references. Either US (8.5x11 inch) or A4 size paper 
can be used. 

13.5 Version 7.30, July 1997 

This version fixes a couple of bugs in the previous version. In particular, the JETTYPE 
selection did not work correctly for supersymmetric Higgs bosons, and there was an error in 
the interactive interface for MSSM input. Since these could lead to incorrect results, users 
should replace the old version. We thank Art Kreymer for finding these problems. 
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Since top quarks decay before they have time to hadronize, they are now put directly 
onto the particle list. Top hadrons (tu, tud, etc.) no longer appear, and FORCE should be 
used directly for the top quark, i.e. 

FORCE 

6,11,-12,5/ 

The documentation has been converted to F/Tp^X. Run either FTj^X 2.09 or FTj^X 2e 
three times to resolve all the forward references. Either US (8.5 x 11 inch) or A4 size paper 
can be used. 

13.6 Version 7.29, May 1997 

While the previous version was applicable for large as well as small tan/3, it did contain 
approximations for the 3-body decays g — > tbWi, 7^ — > bbZj,rrZj, and Wi — > ruZj. The 
complete tree-level calculations for three body decays of the gluino, chargino and neutralino, 
with all Yukawa couplings and mixings, have now been included (thanks mainly to M. Drees). 
We have compared our branching ratios with those calculated by A. Bartl and collaborators; 
the agreement is generally good. 

The decay patterns of gluinos, charginos and neutralinos may differ from previous ex- 
pectations if tan (3 is large. In particular, decays into r's and 6's are often enhanced, while 
decays into e's and /i's are reduced. It could be important for experiments to study new 
types of signatures, since the cross sections for conventional signatures may be considerably 
reduced. 

We have also corrected several bugs, including a fairly serious one in the selection of jet 
types for SUSY Higgs. We thank A. Kreymer for pointing this out to us. 

13.7 Version 7.27, January 1997 

The new version contains substantial improvements in the treatment of the Minimal Su- 
persymmetric Standard Model (MSSM) and the SUGRA model. The squarks of the first 
two generations are no longer assumed to be degenerate. The mass splittings and all the 
two-body decay modes are now correctly calculated for large tan/5. While there are still 
some approximations for three-body modes, ISAJET is now usable for the whole range 
1 ^ tan/3 < M t /M b . The most interesting new feature for large tan/5 is that third genera- 
tion modes can be strongly enhanced or even completely dominant. 

To accomodate these changes it was necessary to change the MSSM input parameters. To 
avoid confusion, the MSSM keywords have been renamed MSSM[A-C] instead of MSSM[l-3], 
and the order of the parameters has been changed. See the input section of the manual for 
details. 

Treatment of the MSSM Higgs sector has also been improved. In the renormalization 
group equations the Higgs couplings are frozen at a higher scale, Q = 'm(£l)M(£r). Run- 
ning t, b and r masses evaluated at that scale are used to reproduce the dominant 2-loop 
effects. There is some sensitivity to the choice of Q; our choice seems to give fairly stable 
results over a wide range of parameters and reasonable agreement with other calculations. 
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In particular, the resulting light Higgs masses are significantly lower than those from Ver- 
sion 7.22. 

The default parton distributions have been updated to CTEQ3L. A bug in the PDFLIB 
interface and other minor bugs have been fixed. 

13.8 Version 7.22, July 1996 

The new version fixes errors in b — > Wt and in some t decays and Higgs decays. It also 
contains a new decay table with updated r, c, and b decays, based loosely on the QQ 
decay package from CLEO. The updated decays are less detailed than the full CLEO QQ 
program but an improvement over what existed before. The new decays involve a number of 
additional resonances, including / (980), <n(1260), / 2 (1270), 1^(1270), 1^(1400), K*(U30), 
Xci,2,3, and ip{2S), so users may have to change their interface routines. 
A number of other small bugs have been corrected. 

13.9 Version 7.20, June 1996 

The new version corrects both errors introduced in Version 7.19 and longstanding errors in 
the final state QCD shower algorithm. It also includes the top mass in the cross sections for 
gb — > Wt and gt — > Zt. When the t mass is taken into account, the process gt — > Wb can 
have a pole in the physical region, so it has been removed; see the documentation for more 
discussion. 

Steve Tether recently pointed out to us that the anomalous dimension for the q — > qg 
branching used in the final state QCD branching algorithm was incorrect. In investigating 
this we found an additional error, a missing factor of 1/3 in the g — > qq branching. The 
first error produces a small but non-negligible underestimate of gluon radiation from quarks. 
The second overestimates quark pair production from gluons by about a factor of 3. In 
particular, this means that backgrounds from heavy quarks Q coming from g — > QQ have 
been overestimated. 

The new version also allows the user to set arbitrary masses for the U(l) and SU(2) 
gaugino mases in the MSSM rather than deriving these from the gluino mass using grand 
unification. This could be useful in studying one of the SUSY interpretations of a CDF 
ee77-^T event recently suggested by Ambrosanio, Kane, Kribs, Martin and Mrenna. Note, 
however, that radiative decay are not included, although the user can force them and multiply 
by the appropriate branching ratios calculated by Haber and Wyler, Nucl. Phys. B323, 267 
(1989). No explicit provision for the decay Z\ — > G7 of the lightest zino into a gravitino or 
goldstino and a photon has been made, but forcing the decay Z\ — > z/7 has the same effect 
for any collider detector. 

A number of other minor bugs have also been corrected. 

13.10 Version 7.16, October 1995 

The new version includes e + e~ cross sections for both SUSY and Standard Model particles 
with polarized beams. The e~ and e + polarizations are specified with a new keyword EPOL. 
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Polarization appears to be quite useful in studying SUSY particles at an e + e" collider. 

The new release also includes some bug fixes for pp reactions, so you should upgrade even 
if you do not plan to use the polarized e + e~ cross sections. 

13.11 Version 7.13, September 1994 

Version 7.13 of ISAJET fixes a bug that we introduced in the recently released 7.11 and 
another bug in g — > qq. We felt it was essential to fix these bugs despite the proliferation of 
versions. 

The new version includes the cross sections for the e + e~ production of squarks, sleptons, 
gauginos, and Higgs bosons in Minimal Supersymmetric Standard Model (MSSM) or the 
minimal supergravity (SUGRA) model, including the effects of cascade decays. To generate 
such events, select the E+E- reaction type and either SUGRA or MSSM, e.g., 

SAMPLE E+E- JOB 
300. ,50000,10,100/ 
E+E- 
SUGRA 

100,100,0,2,-1/ 

TMASS 

170,-1,1/ 

END 

STOP 

The effects of spin correlations in the production and decay, e.g., in e + e~ — > Wi~W[~ , are 
not included. 

It should be noted that the Standard Model e + e~ generator in ISAJET does not include 
Bhabba scattering or W + W~ and Z°Z° production. Also, its hadronization model is cruder 
than that available in some other generators. 

13.12 Version 7.11, September 1994 

The new version includes the cross sections for the e + e~ production of squarks, sleptons, 
gauginos, and Higgs bosons in Minimal Supersymmetric Standard Model (MSSM) or the 
minimal supergravity (SUGRA) model including the effects of cascade decays. To generate 
such events, select the E+E- reaction type and either SUGRA or MSSM, e.g., 

SAMPLE E+E- JOB 
300. ,50000,10,100/ 
E+E- 
SUGRA 

100,100,0,2,-1/ 

TMASS 

170,-1,1/ 

END 

STOP 



81 



The effects of spin correlations in the production and decay, e.g., in e + e~ — > W{~ , are 
not included. 

It should be noted that the Standard Model e + e~ generator in ISAJET does not include 
Bhabba scattering or W + W~ and Z°Z° production. Also, its hadronization model is cruder 
than that available in some other generators. 

13.13 Version 7.10, July 1994 

This version adds a new option that solves the renormalization group equations to calculate 
the Minimal Supersymmetric Standard Model (MSSM) parameters in the minimal super- 
gravity (SUGRA) model, assuming only that the low energy theory has the minimal particle 
content, that electroweak symmetry is radiatively broken, and that R-parity is conserved. 
The minimal SUGRA model has just four parameters, which are taken to be the common 
scalar mass mo, the common gaugino mass mi/2, the common trilinear SUSY breaking term 
Aq, all defined at the GUT scale, and tan (3; the sign of [i must also be given. The renormal- 
ization group equations are solved iteratively using Runge-Kutta integration including the 
correct thresholds. This program can be used either alone or as part of the event generator. 
In the latter case, the parameters are specified using 

SUGRA 

m , mi/2, A , tan/5, sgn/i 

While the SUGRA option is less general than the MSSM, it is theoretically attractive and 
provides a much more managable parameter space. 

In addition there have been a number of improvements and bug fixes. An occasional 
infinite loop in the minimum bias generator has been fixed. A few SUSY cross sections 
and decay modes and the JETTYPE flags for SUSY particles have been corrected. The 
treatment of B baryons has been improved somewhat. 
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